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16.  Abstract 
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SUMMARY 


This  report  presents  the  results  of  a  program  of  analysis  and  test 
to  investigate  the  effects  of  temperature,  dissipation  level  and 
geometry  on  the  dynamic  behavior  of  elastomer  elements.  The  report 
presents,  first,  a  review  of  force  displacement  relationships  in 
elastomer  elements  and  the  effects  of  frequency,  geometry  and 
temperature  upon  these  relationships.  Based  on  this  reviev;,  methods 
of  reducing  stiffness  and  damping  data  for  shear  and  compression  test 
elements  to  material  properties  (storage  and  loss  moduli)  and 
empirical  geometry  factors  are  developed  and  tested  using  previously 
generated  experimental  data.  A  prediction  method  which  accounts  for 
large  amplitudes  of  deformation  is  developed  on  the  assumption  that, 
their  effect  is  to  increase  temperature  through  the  elastomers,  tlie^^e- 
by  modifying  the  local  material  properties.  Various  simple  methods 
of  predicting  the  radial  stiffness  of  ring  cartridge  elements  are 
developed  and  compared. 

The  test  method  for  all  dynamic  stiffness  and  damping  measurements 
involves  resonating  a  mass,  supported  on  the  test  element,  by  exciting 
the  base  to  which  the  test  element  is  mounted.  The  specimens  tested 
were:  Two  rectangular  shear  specimens;  one  cylindrical  compression 
specimen;  and  two  cartridge  elements  of  different  length.  Tests  were 
performed  at  temperatures  between  5  and  93  C  and  at  dissipation  levels 
between  0.022  watts/cm^  and  4.4  watts/cm^  material. 


Material  properties  were  determined  from  the  shear  specimen  tests  as 
a  function  of  frequency  and  temperature.  Using  these  material 
properties,  numerical  predictions  of  stiffness  and  damping  for 
cartridge  and  compression  specimens  were  made  and  compared  with  corre¬ 
sponding  measurements  at  different  temperatures,  with  encouraging 
results . 

Successful  comparisons  of  high  dissipation  predictions  and  measure¬ 
ments  were  made  for  the  compression  specimens  but  it  is  shown  that, 
for  the  shear  specimen,  strain  is  an  important  influential  parameter, 
even  when  the  temperature  rise  is  small. 

The  method  of  reduced  variables  (which  allows  frequency  dependent 
data  at  different  temperatures  to  be  collapsed  to  a  single  curve  in 
reduced  variables)  is  tested  out  on  the  material  property  results. 

It  is  shown  that  the  storage  modulus  follows  the  reduced  variable  law 
satisfactorily  whereas  the  loss  modulus  does  not. 


INTRODUCTION 


Tlic  use  of  support  damping  as  a  means  to  control  rotating  machinery 
vibrations  is  seeing  increasingly  wide  application  in  advanced  turbo¬ 
machinery.  It  will  also  have  an  important  role  to  play  in  advanced, 
flexible,  power  transmission  shafting.  Presently  the  most  common 
type  of  damper  in  these  applications  is  the  squeeze-film  damper  in 
parallel  with  some  type  of  mechanical  flexure. 

Elastomer  dampers  are  an  attractive  alternative  to  the  squeeze  film 
for  rotating  machinery  and  other  applications  because  of  their  sim¬ 
plicity^  their  inherent  combination  of  stiffness  and  damping*  their 
compactness;  and  their  lack  of  need  for  seals  or  oil  supply.  In  the 
form  of  0-rings  or  cartridges  they  are  being  considered  for  low-cost 
engine  applications  and  for  helicopter  transmission  shafting. 

Two  of  the  factors  which  resist  the  growth  in  application  of  elastomer 
dampers  are  the  limited  availability  of  design-oriented  data  on  their 
dynamic  behavior  and  limited  quantification  of  the  problems  to  be 
encountered  in  their  application.  Dynamic  testing  under  controlled 
conditions,  coupled  with  practical  interpretation  of  the  test  results 
is  the  means  to  fill  this  need.  Thereby,  the  influence  of  important 
geometrical,  environmental  and  chemical  design  parameters  can  be 
determined , 

A  program  for  testing  and  analysis  of  elastomer  materials  has  been 
in  progress  at  MTI  since  19/1.  The  objectives  of  the  program  are  to 
develop  elastomer  dynamics  Leclmology  and  to  reduce  the  technology 
to  a  form  in  which  it  can  be  readily  applied  in  design. 

Experimental  and  analytical  methods  are  employed  in  meeting  these 
objectives.  Firstly,  from  a  limited  number  of  tests,  general  prop¬ 
erties  are  sought  which  characterize  an  elastomer  material.  Secondly, 
prediction  methods  are  sought  which  employ  these  material  properties, 
in  combination  with  appropriate  geometry  effects,  to  yield  values 
for  stiffness  and  damping  of  particular  elastomer  configurations. 

This  report  is  the  third  in  a  series.  Preceding  reports  (refs.  1 
and  2)  have  described  the  development  and  refinement  of  a  powerful 
test  technique,  the  Base  Excitation  Resonant  Mass  (BERM)  method, 
and  its  initial  application  to  the  development  of  elastomer  element 
properties  in  a  controlled  temperature  environment  under  conditions 
of  low  dissipation.  The  BERM  method  most  conveniently  supplies  data 
in  the  100  Hz  to  1000  Hz  frequency  range  -  into  which  the  majority 
of  rotating  machinery  vibration  problems  fall.  Important  features 
of  the  method  are  its  controlled  temperature  environment;  its  ability 
to  apply  controlled  high  or  low  levels  of  dissipation  at  frequencies 
up  to  1000  Hz  or  higher;  its  exploitation  of  mechanical  resonance  to 
achieve  phase  angles  clustered  around  90  degrees;  and  its  determination 


oi'  stiffness  and  damping  from  accelerometer  signals  only,  each 
processed  by  identical  electronics,  to  minimi2:e  spurious  phase 
errors  in  the  data. 

The  work  presented  herein  provides  a  much  expanded  body  of  test 
data  under  a  range  of  frequencies,  temperatures  and  dissipation 
levels;  development  of  interpretive  and  predictive  methods;  appli^ 
cation  of  these  methods;  and  verification  of  their  ability  to 
predict  dynamic  behavior.  As  a  framework  against  which  this  work 
can  be  viewed,  Figure  1  is  presented. 

Figure  1  is  a  statement  of  the  possible  interactions  between  tesc3 
data  reduction,  analysis,  performance  prediction  and  verification. 

As  shown,  the  starting  point  for  practical  information  is  testing, 
and  the  BERM  test  method  provides  the  main  source  of  test  data 
within  this  program.  The  output  of  the  tests  is  raw  data  (ampli¬ 
tudes,  phase  angles,  etc)  which  is  converted  by  automated  data 
analysis  methods  to  specimen  properties  -  that  is,  stiffness  and 
damping  of  the  elements  being  tested,  as  a  function  of  frequency. 

To  replace  the  discrete  data  points  resulting  from  tests  with  a 
continuous  relationship  which  can  be  handled  mathematically,  an 
appropriate  frequency  correlation,  preferably  of  minimum  complexity, 
is  performed.  The  output  of  this  correlation  provides  a  compact 
definition  of  the  dynamic  behavior  of  the  elements  being  tested 
which,  together  with  similar  compact  definitions  for  other  elements, 
yields  a  specimen  catalog;  that  is,  a  list  of  stiffness  and  damping  . 
properties  of  a  series  of  different  elements  under  different  con¬ 
ditions  of  temperature  and  power  level. 

To  provide  more  universally  applicable  information,  the  specimen 
information  is  processed  to  yield  material  properties  (which  relate 
stress  to  strain)  as  a  function  of  frequency.  Where  available, 
other  sources  of  test  data  may  also  be  employed.  The  output  of  this 
step  is  a  material  properties  catalog,  describing  one  material  or 
several  materials  as  a  function  of  frequency  and  temperature. 

To  apply  elastomer  material  properties  (moduli)  in  the  design  pre¬ 
diction  of  elastomer  element  properties,  it  is  necessary  to  trans¬ 
late  stress-strain  relationships  back  into  force  -  displacement 
relationships.  For  each  elastomer  configuration,  there  exists  a 
ratio  between  the  stiffness  and  damping  at  a  particular  frequency 
and  temperature  and  the  moduli  of  the  material  at  the  same  frequency 
and  temperature.  These  ratios,  or  geometry  factors,  will,  in 
general,  be  a  function  of  frequency,  temperature,  amplitude  and 
geometry . 

As  shown  in  Figure  1,  sources  of  these  ratios  are  solutions  to  the 
continuum  equations  for  the  elastomer  element,  in  some  cases  com¬ 
bined  with  empirical  coefficients  obtained  by  tests.  The  spectrum 
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of  solution  complexity  can  run  from  a  simple  algebraic  formula  to 
the  solution  of  complex  continuum  energy  and  viscoelasticity  equa¬ 
tions  in  two  or  three  dimensions.  The  body  of  ratios,  or  geometry 
factors,  determined  in  one  or  several  ways,  is  the  geometry  catalog. 

The  material  properties  catalog  and  the  geometry  catalog  are  the  two 
sources  of  information  which  embody  the  predictive  model.  They  make 
it  possible  to  predict  stiffness  and  damping  as  a  general  function 
of  material,  geometry,  frequency,  amplitude  and  temperature. 

Progri'ss  along  most  of  the  paths  of  Figure  1  has  been  made  and  r 
this  program,  as  will  be  described  in  this  report. 

Important  questions  to  be  addressed  by  this  technology  program  are. 

"What  level  of  complexity  and  sophistication  is  needed  to 
generate  an  adequate  geometry  catalog  for  design  purposes?" 

"what  are  the  benefits,  in  improved  accuracy,  to  be  expected 
from  advanced  solutions  relative  to  simple  solution  methods?" 

The  simplest  methods  of  all  provide  geometry  factors  which  are  in¬ 
dependent  of  frequency,  and  in  two  cases  such  predictive  methods 
have  been  tested  and  their  accuracy  assessed  against  test  data  or 
more  complete  representations.  These  methods  clearly  lack  rigor  and 
there  are  discrepancies  between  their  predictions  and  measured  data. 

At  the  same  time,  for  certain  application  of  elastomers,  their  pre¬ 
cision  may  be  adequate.  In  other  applications,  the  development  of, 
as  yet,  unavailable  sophisticated  prediction  methods  will  be  required. 
Given  the  level  of  component  prediction  accuracy  to  be  expected,  as 
indicated  in  this  report,  sensitivity  studies  for  particular  applica¬ 
tions  will  reveal  the  likely  uncertainties  in  predicted  system  per¬ 
formance  which  will  result. 

The  following  sections  of  the  report  present,  first,  a  summary  of 
results,  conclusions,  and  recommendations,  including  the  material 
properties  and  geometry  catalog  resulting  from  the  tests  and  anal¬ 
ysis  presented  in  the  report.  Then  follows  a  review  of  force- 
displacement  relationships  for  elastomers;  the  development  of 
interpretive  and  predictive  models;  a  description  of  test  methods, 
including  parameter  ranges;  a  presentation  of  the  test  results;  and 
a  discussion  of  these  results. 
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SUMMARY  OF  RESULTS 


'i 

s 


1.  Methods  of  identifying  generalized  material  properties  and 
their  relationship  to  dynamic  behavior  of  elastomer  elements 
have  been  developed,  applied  and  tested. 

2.  Dynamic  test  data  for  several  geometries  has  been  generated, 
for  dissipation  levels  from  0.022  to  4,4  watts/cm^j  for  tem¬ 
peratures  from  5  to  93  C,  and  for  frequencies  from  70  to 
1500  Hz . 

3.  The  geometries  considered  are; 

•  shear  loaded  rectangular  elastomer  strips, 

•  compression  loaded  right  solid  cylinders,  and 

•  radially  loaded  cylindrical  cartridges. 

4.  The  material  is  a  broad  temperature  range  (BTR)  elastomer - 
polybutadiene . 

5.  Material  property  values  (storage  and  loss  moduli)  have  been 
generated,  as  a  function  of  frequency  and  temperature,  from 
the  shear  elements. 

6.  Geometry  factor  values  (analogous  to  the  traditional  shape 
factor  coefficient)  have  been  generated  for  compression  loaded 
right  circular  cylinders  as  a  function  of  frequency. 

7.  These  geometry  factors  have  been  shown  to  provide  a  self- 
consistent  representation  of  data  generated  under  previously 
reported  work. 

8.  These  geometry  factors  have  provided  reasonably  effective 
predictions  of  stiffness  and  damping  for  a  range  of  temperatures 
other  than  those  for  which  they  were  developed. 

9.  These  empirical  geometry  factors,  which  depend  on  frequency, 
have  been  compared,  in  terms  of  their  ability  to  predict 
measured  stiffness  and  damping,  with  use  of  a  frequency  in¬ 
dependent  static  shape  factor  coefficient. 

10.  Cartridge  test  specimens  and  an  appropriate  test  fixture  have 
been  designed  and  fabricated. 

11.  Predictions  of  cartridge  stiffness  and  damping  have  been  made 
as  a  function  of  frequency,  using  two  different  algebraic 
expressions.  The  simplest  of  these  is  based  on  representation 
of  the  cartridges  as  a  series  of  noninteracting  beams  and 
column  elements, 

12.  These  two  simple  methods  predict  results  which  differ  by  a 
factor  of  two  but  provide  a  band  within  which  the  great 
majority  of  measured  values  lie. 
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13.  A  thermoviscoelastic  prediction  method  indicates  some  similar 
trends  for  stiffness  as  observed  under  conditions  of  high  dis-* 
sipation;  that  is,  decreasing  stiffness  with  increasing  dis¬ 
sipation  levels.  This  statement  applies  particulary  to  compres¬ 
sion  test  specimens. 

14.  Although  the  conductivity  is  unknown  for  this  material,  an  effective 
thermal  conductivity  value  has  been  generated  which  gives  a 

good,  quantitative  match  between  measured  and  predicted  values  of 
high  dissipation  stiffness  and  damping  for  compres^^ion  loaded  cylinders. 

15.  Use  of  the  same  effective  conductivity  value  for  shear  eleii:'its 
gives  predicted  stiffness  and  damping  values  which  are  higher 
than  observed,  under  high  dissipation  conditions. 

16.  Measured  stiffness  and  damping  values  were  correlated  against 
strain  and  against  dissipation.  More  consistent  dependence 
on  strain  than  on  dissipation  was  observed. 

17.  A  "strain  softening"  effect  is  apparent  in  the  shear  specimen 
results . 


18.  The  method  of  reduced  variables  was  effective  in  reducing 
storage  modulus  data  for  different  temperatures  to  a  single 
curve  of  reduced  modulus  vs  reduced  frequency,  but  was  not  as 
effective  for  loss  modulus.  Strain  dependence  of  the  results 
is  shown  to  be  a  likely  cause  of  this  apparent  inconsistency 
for  loss  modulus . 

19.  The  material  property  catalog  for  the  tested  elastomer,  based 
on  the  results  presented  in  this  report  is  as  follows: 


Temper- 

ature 

32  C 
43  C 
66  C 
93  C 


1  Storage  Modulus 

G^(aj) 

=  3.686 

X 

10^ 

0.2037 

0) 

=  1.904 

X 

10^ 

0.2627 

u> 

G'M 

=  1.097 

X 

10^ 

0.3037 

CO 

G'Oi.) 

=  1.197 

X 

10^ 

0.2450 

cu 

Loss  Modulus 


G"(u))  =  8.333 

X 

10^ 

-0.1277 

0) 

G"(ui)  =  3.971 

X 

10^ 

-0.0637 

0) 

G"tu)  =  2.423 

X 

10^ 

-0.0301 

(1) 

G^(a))  =  1.064 

X 

o 

o^ 

+0.2246 

U) 

20.  The  geometry  catalog  for  the  three  tested  goemetries  is  as 
follows : 
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CONCLUSIONS  AND  RECOMMENDATIONS 


On  the  basis  of  the  results  presented  in  this  report,  the  following 

conclusions  are  drawn. 

1.  Storage  and  loss  modulus  values,  obtained  as  a  function  of 
frequency  from  shear  specimen  tests,  are  an  effective 
representation  of  the  dynamic  properties  of  a  material. 

2.  Empirical,  frequency  dependent  coefficients  provide  a  self- 
consistent  representation  of  the  variation  of  stiffness  ith 
shape  factor  for  cylindrical  compression  elements. 

3.  For  new  conditions,  the  frequency  dependent  empirical  co¬ 
efficients  provide  similar  predictive  ability  to  the  use  of  a 
constant,  static,  coefficient. 

4.  Either  method  of  prediction  for  the  compression  elements 
incurs  errors  of  10  to  30  percent  (high)  relative  to  test 
data  generated  under  the  present  program. 

5.  Predictive  ability,  based  on  simple  static  analyses,  for 
radially  loaded  ring  cartridges  provides  upper  and  lower  limits 
for  stiffness  and  damping  within  which  the  great  majority  of 
measured  data  lies. 

6.  The  effect  of  high  dissipation  is  to  reduce  stiffness. 

7 .  According  to  the  frequency  and  geometry  damping  can  be  re¬ 
duced  or  increased  by  high  dissipation. 

8.  Self-heating  and  strain-softening  effects  are  the  combined 
causes  of  stiffness  and  damping  changes  under  high  dissipa¬ 
tion. 

9.  For  all  results  presented,  stiffness  and  damping  depend  more 
consistently  on  strain  than  on  dissipation  rate. 

10.  Increasing  strain  up  to  0.02  (double  amplitude)  generally  in¬ 
creases  the  loss  coefficient  for  elements  of  polybutadiene. 

For  future  investigations,  the  following  recommendations  are  made. 

1.  The  understanding  of  the  dynamic  behavior  of  elastomer  elements 
will  be  improved  by  obtaining  further  data  at  very  low  strain 
levels,  and  by  relating  it  to  the  present  data  at  high  strains. 

2.  Simple  methods  for  including  strain  dependence  of  dynamic  be¬ 
havior  in  elastomer  design  practice  should  be  developed. 


3*  The  thermal  properties  of  the  elastomer  material  under  investi¬ 
gation  should  be  determined. 

4.  Material  and  geometrical  property  tests  for  a  number  of  dif¬ 
ferent  elastomer  materials  should  be  performed  to  determine  the 
extent  to  which  the  presently  reported  results  can  be 
generalized. 

5.  Predictive  ability  for  radially  loaded  cartridges  must  be 
further  evaluated  in  a  rotating  load  environment. 

6.  Because  they  are  compact,  readily  available,  elastomer 
elements,  the  dynamic  properties  of  0-rings  should  be  investi¬ 
gated  by  analysis  and  by  test,  with  an  effort  being  made  to 
relate  the  dynamic  properties  to  the  basic  material  properties 
of  storage  and  loss  moduli. 

7.  Application  studies,  analytical  then  experimental,  should  be 
performed  to  determine  the  benefits  and  limitations  of 
elastomer  elements  as  vibration  control  elements. 

8.  Successful  short-term  application  tests  should  be  followed  by 
tests  for  life  and  survivability  to  determine  the  extent  to 
which  these  problems  will  limit  application. 

9.  Application  studies  should  be  interpreted  to  reveal  what  level 
of  predictive  accuracy  is  required  for  elastomer  elements. 

10.  More  sophisticated  prediction  methods  should  be  developed  as 
indicated  by  the  application  studies. 
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REVIEW  OF  MATERIAL  PROPERTIES  AND  THE  EFFECTS 
OF  GEOMETRY,  FREQUENCY,  AMPLITUDE  AND  TEMPERATURE 


In  this  section  the  relationships  between  force,  displacement,  stress 
and  strain  in  an  elastomer  specimen  are  reviewed,  with  particular 
emphasis  on  frequency  domain  behavior  and  the  influence  of  geometry, 
temperature  and  amplitude. 

Relaxation  and  Hysteresis 

Elastomer  deformation  properties  are  dominated  by  the  process  of  re¬ 
laxation.  By  this  process,  the  deflection  for  a  given  load  (strain 
for  a  given  stress)  increases  as  a  function  of  time  »  tendin>  towards 
an  asymptotic  limit  after  a  period  of  sufficient  length.  For 
short  periods  of  time,  the  displacement  under  a  particular  load  wi.i 
be  significantly  lower  than  the  sympototic  limit. 

If  a  load  cycle  is  applied  repeatedly  to  an  elastomer  element,  the 
process  of  relaxation  causes  hysteresis  -  that  is,  a  load  deformation 
curve  which  depends  on  the  rate  of  change  of  load  as  well  as  on  the 
load's  magnitude  -  see  Figure  2,  in  which  a  representative  elliptical 
hysteresis  loop  is  presented.  Implicit  in  such  a  load  deflection 
curve  or  hysteresis  loop  is  a  loss  of  energy  per  cycle  equal  to  the 
area  of  the  hysteresis  loop.  Expressed  in  vibration  terminology, 
elements  with  such  a  load-deflection  curve  provide  dissipation,  or 
damping.  As  shown  by  Lazan  (ref.  3),  for  an  elliptical  hysteresis 
loop,  the  material  behavior  is  linear,  and  the  damping  may  be  expressed 
as 


D 

K^Ccd)  =  (1) 


where  D^  is  the  area  of  the  hysteresis  loop  in  energy  units  and  Y  is 
the  single  amplitude.  It  is  also  possible  to  define  stiffness  in^ 
terms  of  this  loop  as  the  ratio  between  force  at  maximum  delfection 
and  the  value  Y^  of  maximum  deflection. 

=  F^/Y^  (2) 

where  F^  and  Y^^  arc  defined  in  Figure  2.  In  the  following  paragraphs 
the  relationships  between  instantaneous  force  and  displacement  are 
developed . 

It  is  noted,  at  this  point,  that  the  linearity  of  elastomer  be- 
liavior  has  been  investigated  by  Warnaka  (ref.  4)  for  materials 
including  neoprene,  nitride  rubber  and  BTR  elastomer  for  strains 
up  to  0.1.  Using  electronics  to  sense  deviations  from  an  el¬ 
liptical  hysteresis  loop,  Warnaka  concluded  that  linear  dynamic 


behavior  within  a  small  discrepancy  was  exhibited  by  elastomer  elements 
for  even  the  highest  strain  at  frequencies  in  the  100-1000  Hz  range. 

Relationships  Between  Force  And  Displacement  in  The  Frequency  Domain 

A  linear  elastomer  element  undergoing  a  repeated  deflection  cycle 
defined  by 


Y 


V* 

Y  e 


(3) 


gives  rise  to  a  cyclic  force  variation  described  by 


F 


i^jDt 

F  e 


(4) 


where: 


O) 

Y 

F 

*  h 

Y  ,F 


is  the  frequency  (rad/sec) 
is  the  instantaneous  displacement  (m) 
is  the  instantaneous  force  (N) 
is  time 

are  complex  numbers  defining  the  amplitude  and  phase 
of  the  displacement  and  force  variations,  (and  it  is 
understood  that  only  the  real  parts  of  expressions 
(3)  and  (4)  are  used  to  establish  values  for  Y  and  F). 


The  force  and  displacement  may  be  related  by 
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F  =  K  Y 


where  K 


* 


is 


a  complex  number 


* 

K 


+  iK^ 


(5) 

(6) 


and  is,  generally,  a  function  of  frequency.  The  quantities  and 

will  be  referred  to  as  the  stiffness  and  dmaping  of  the  elastomer 

element.  The  K-  and  K^,  which  depends  on  frequency,  temperature 

amd  amplitude  characterize  the  dynamic  behavior  of  the  element.  The 

damping  ability  of  an  elastomer  element  is  sometimes  expressed  in 

terms  of  a  loss  coefficient  (nondimensional)  which  relates  damping 

to  stiffness,  or  the  energy  dissipated  percycle  (D  )  to  the  energy 

stored  at  maximum  displacement  (U  ) 

s 


■^2  <“>  _ 

Kj^  (od)  2  n  Ug 


(7) 


Relationships  Between  Stress  And  Strain 

Considering  any  elastomer  element,  its  geometry  can,  normally,  be 
characterized  by  a  certain  area  of  the  element  (A)  which  is  stressed 
and  a  certain  dimension  of  the  element,  h,  which  is  strained. 

Figure  3  shows  two  elastomer  element  geometries  which  have  been 


LOAD 


WIDTH, 

-  b  - 


)  1  i  I 
MM 
MM 
MM 
U-X-i 


LENGTH, 


^thickness,  h 


A.  SHEAR  ELEMENT 


LOAD 


B  COMPRESSION  ELEMENT 


Tost  [‘llcmonts  l'’or  Shear  and  Compression  Loading 


tested  for  dynamic  stiffness  and  dampiiig,  and  which  illustrate  the 
effects  of  geometry.  For  the  shear  element  (Fig.  3A) ,  the  stressed 
area,  A,  is  its  width  times  its  length  (b-x)  and  the  strained  dimen¬ 
sion  is  its  tliickness,  h.  For  the  compression  element  (Fig.  3B), 

I  he  stressed  area  is  the  top  area  of  the  cylinder  (tt  D^/4)  and  its 
strained  dimension  is  the  height,  h. 

For  the  shear  element  stiffness  and  damping  can  be  related  to  effec¬ 
tive  moduli  as  follows. 


Stiffness : 

■'i  ■  Kfct 

(8) 

Damping: 

"2  ■  “^eff  h 

where  are  shear  moduli,  commonly  referred  to  as  the 

effective  storage  and  loss  moduli,  respectively. 


For 

the 

the  compression  element  stiffness 
effective  moduli  as  follows. 

and  damping 

can  be  related 

Stiffness ; 

‘'i  = 

3  G'  ..  - 

eft  h 

(10) 

Damp ing : 

K 

2 

3  g"  ^ 

eff  h 

(11) 

where  the  factor  of  3  arises  from  the  classical  elasticity  relation¬ 
ship  E  =  2(1  +  v)G  and  the  fact  that  elastomer  Poisson's  ratios  are 
typically  in  the  range  0.4945  to  0.4999  (Holownia,  ref.  5).  Any 
errors  in  applying  this  classical  elasticity  relationship  to  dynamic 
behavior  of  a  viscoelastic  material  will  be  implicit  in  the  state¬ 
ment  that  ate  effective  moduli.  G  g" 

functions  of  the  ® 

the  temperature,  the  initial  strain  and  the  geometry.  They  are, 
therefore,  component  rather  than  material  properties.  Further  dis¬ 
cussion  will  explore  the  relationships  between  effective  and  "true" 
moduli , 


are,  generally 
material,  the  frequency,  the  amplitude  of  loading, 


The  Effects  of  Geometry 

For  static  loading  the  effects  of  bending  can  be  shown  (ref.  6)  to 
cause  the  following  relationship  between  effective  and  true  modulus 
for  the  shear  element  of  Figure  3A. 


<":ff>o 


=  G 


_1 


+  h2/3x2  J 


(12) 


where  X  is  the  dimension  of  the  sheared  area  in  the  direction  parallel 
to  the  loading,  h  is  the  thickness,  and  the  subscript,  o,  implies 
static  loading.  In  Figure  4,  the  stiffness  correction  factor 


15 


STIFFNESS  CORRECTION  FACTOR 


MAXIMUM  (h/x)  FOR  — 
LESS  THAN  1%  COMBINED 
STRESS  EFFECTS 


0  2  .4  .6  .0  1.0  I 

(h/x)-  RATIO  OF  THICKNESS  TO  LENGTH  PARALLEL  TO  DIREC 
TION  OF  LOADING 


Static  Stitfncss  Correction  Factor 
Spec Imen 


Reel angu lar  Shear 


is  plotted  as  a  function  of  X/h  for  a  rectangular  shear 
clement.  Under  static  compression  loading,  the  interaction  of  axial 
radial  and  circumferential  stresses  was  shown  by  Gent  and  Lindley 
(ref.  7),  Hatiorl  and  Takei  (ref.  8),  and  Payne  (ref.  9)  to  give  the 
following  relationship  between  apparent,  or  effective  static  modulus 
and  actual,  or  true,  static  modulus: 

'“e'ff'o  -  (13) 


where  S  is  a  shape  factor  defined  as  the  ratio  of  the  loaded,  .area 
to  the  unloaded  area  (D/4h  for  a  cylinder)  and  8  is  a  coefficx  ;Pt, 
shown  by  Gent  and  Lindley  (ref.  7)  to  be  equal  to  2 . 

A  more  complete  expression  based  on  the  work  of  Gent  and  Lindley 
(ref.  7),  which  attempts  to  account  for  bulk  compressibility  effects 
is : 


G  3  0' 

7g )  =  (1  +  (14) 

^  eff^o  ^ 

where  B  is  the  bulk  modulus  for  the  material.  However,  this  ex¬ 
pression  has  been  tested  by  Moghe  and  Neff  (ref.  10)  using  classical 
elasticity  solutions  and  shown  also  to  be  incomplete  for  large 
values  of  S  in  that  an  additional  dependence  on  Poisson's  ratio 
occurs.  For  values  of  S  less  than  10,  however,  neither  the  second 
term  in  Eq .  (14)  nor  the  inadequacies  due  to  Poisson's  ratio  effects 
should  be  significant. 

In  Figure  5,  the  static  stiffness  correction  factor  ( (<^ef f ^o/^o\ 
as  given  by  Eq .  (13)  with  p  =  2,  is  plotted  against  the  ratio  of 
diameter  to  height,  for  a  compression  loaded  cylinder 

The  existence  of  such  a  square  law  dependence  on  shape  factor  has 
been  verified  experimentally  by  Payne  (ref.  9).  He  showed  that,  for 
natural  rubbers,  the  shape  coefficient,  p,  was  close  to  2  but,  for 
harder  materials  the  coefficient  p  could  fall  by  as  much  as  50  per¬ 
cent  , 

Under  low  frequency  dynamic  loading  (10  Hz)  Payne  also  demonstrated 
that  the  effective  storage  modulus  varied  with  shape  factor  according 
to  a  square  law. 

The  analytical  or  experimental  quantification  of  the  effects  of 
geometry  (shape)  at  higher  frequencies  appears  to  be  limited.  Nashif 
Cannon  and  Jones  (ref.  11)  have  demonstrated  that  shape  has  a  pro¬ 
nounced  influence  on  the  effective  storage  modulus  of  cylindrical 
specimens,  in  the  frequency  range  of  250  to  500  Hz,  but  have  not 
generalized  their  results. 
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The  Effects  of  Frequency 

As  might  be  expected  from  a  material  whose  behavior  is  governed  by 
time  dependent  relaxation  processes,  the  slope  and  area  of  the 
hysteresis  loop,  together  with  the  stiffness,  damping,  and  loss 
coefficient  are  functions  of  the  frequency  with  which  the  load 
cycle  is  executed.  It  is  generally  observed  that  the  dynamic 
stiffness  is  higher  than  the  static  value,  and  that  stiffness  tends 
to  increase  with  increasing  frequency.  Shown  on  a  scale  which  may 
cover  many  decades  of  frequency,  the  general  form  of  elastomer  be¬ 
havior  variation  has  three  regions  (ref.  12): 

•  Rubbery 

•  Transition  (Leathery) 

•  Glassy 

These  are  illustrated  in  Figure  6, 

In  the  rubbery  region,  the  behavior  is,  as  implied  **rubbery.”  Elastomer 
elements  are  resilient  and  the  damping  is  usually  of  the  order  of  10 
to  30  percent  of  the  stiffness.  The  hysteresis  loop  is  a  rather 
narrow  ellipse.  As  the  frequency  increases,  the  transition  region 
is  encountered.  Both  stiffness  and  damping  increase  more  rapidly  in 
this  region  but  the  damping  moves  ahead  of  the  stiffness  and  the 
hysteresis  ellipse  becomes  broader. 

For  a  limited  frequency  range,  damping  (K2)  becomes  higher  than 
stiffness  (K]^)  and  the  loss  coefficient  becomes  greater  than  1.  For 
further  increases  in  frequency  the  stiffness  increases  more  than  the 
damping  and  in  the  glassy  region  both  quantities  level  off  at  values 
giving  a  loss  factor  which  again  tends  to  be  0.3  or  less. 

The  rubbery  and  transition  regions  are  those  normally  encountered  in 
elastomers  and  offer  the  best  opportunity  for  use  in  vibration  con¬ 
trol  . 

Representation  of  Frequency  Dependence 

Observed  variation  of  stiffness  and  damping  with  frequency  gives 
rise  to  discrete  data  points.  For  mathematical  purposes  it  is  con¬ 
venient  to  seek,  for  an  element,  a  continuous  relationship  between 
force  and  displacement  (or  stress  and  strain)  which  closely  matches 
observed  behavior.  There  are  several  possible  approaches  to  this 
problem,  with  different  levels  of  sophistication  and  convenience, 
as  discussed  below. 

1_.  Generalized  Viscoelastic  Models.  An  approach  to  representing 
the  variation  of  stiffness  and  damping  of  an  elastomer  element  with 
frequency  is  to  replace  the  elastomer  with  a  mechanical  analog,  con¬ 
sisting  of  a  series  of  springs  and  viscous  dashpots.  The  simplest 
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of  these  models  are  the  Maxwell  and  Voigt  elements  (Fig,  7).  Using 
these  basic  elements  as  building  blocks,  more  complex  systems  of 
springs  and  dashpots  may  be  developed,  such  as  those  shown  in 
Figure  8, 

The  application  of  a  triparameter  model,  as  shown  in  Figure  9,  to 
elastomer  specimen  test  data  generated  by  the  BERM  method  was  dem¬ 
onstrated  by  Chiang,  Tessarzik  and  Badgley  (ref.  1).  The  methodology 
to  determine,  for  a  selected  general  viscoelastic  model,  the  element 
values  to  best  fit  a  particular  variation  of  stiffness  and  damping 
with  frequency  was  developed  and  demonstrated  by  Gupta,  Tes.'^rzik 
and  Cziglenyi  (ref.  2), 

2,  Generalized  Force  -  Displacement  Relationships.  An  a  Iternat ve 
approach  for  representing  the  variation  with  frequency  is  to  assume 
a  general  linear  differential  relationship  between  force  and  dis¬ 
placement: 


2  3 

_  .  dF  .  d  F  ,  d^^F  , 

V  “l  dT  +  “2  77  +  "3  73 

dt  dt 


=  c  Y  + 
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dY  ^ 

"1  ^  + 


(15) 


For  sinusoidal  loading  and  displacement: 


F  = 

Y  = 

So  that  Eq .  (15)  becomes: 
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A  close  connection  between  this  generalized  force  displacement  re¬ 
lationship  and  the  generalized  viscoelastic  models  exists  since,  for 


;i  given  order  of  expansion  on  left-  and  right-hand  sides  of  Eq .  (15), 
(here  exists  a  corresponding  combination  of  springs  and  dashpots 
(ref.  i2),  which  will  give  the  same  I orce-disp lacement  relationship. 


For  example  the  differential  equation  for  the  triparameter  model 
(Fig.  9)  is 


Ft  .11  ‘1*’’ 

Si?' 


c. 


K, 


"  (^,.2  "t,*  f 


(20) 


so  that  the  nonzero  coefficients  of  Eq  .  (15)  for  the  tr iparaiu.  t er 
mo(iel  are 


a  =  1 

o 


Cti/Kti 


=  K 


T2 


'I  Si  ^*^12 


(21) 


3.  Direct  Fit  of  Mathematical  Expressions.  A  direct  approach  to 
the  problem  is  to  impose  upon  the  data  some  correlation  function  and 
to  seek  those  coefficients  which  give  the  best  match  between  the 
function  and  the  data.  Polynomial,  trigonometrical,  or  power  law 
(Ki  =  relationships  are  some  of  the  options  available  for  this 

approach.  It  was  found  by  Gupta,  et  al  (ref.  2),  that  the  use  of  a 
power  law  relationship  was  quite  adequate  to  represent  the  behavior 
of  elements  of  polybutadiene  in  the  frequency  range  100  to  1000  Hz. 
For  engineering  purposes,  a  known  relationship  between  stiffness, 
damping  and  frequency  for  elastomer  elements  of  different  geometries 
is  all  that  is  required  to  use  these  elements  in  design  for  control 
of  machinery  vibrations.  In  essence  this  approach  may  be  considered 
as  the  use  of  a  generalized  Voigt  element  in  which  the  storage  and 
dissipative  elements  are  a  continuous  function  of  frequency. 


Influence  of  temperature.  Elastomer  behavior  is,  in  general, 
strongly  influenced  by  temperature.  Thus  complete  characterization 
of  elastomer  elements  must  provide  their  dynamic  behavior  over  the 
frequency  and  temperature  ranges  of  interest.  In  considering  the 
effects  of  temperature  the  following  observations  are  often  made. 


1.  Temperature  tends  to  be  most  influential  where 
frequency  is  also  strongly  Influential. 

2.  Increasing  temperature  tends  to  produce  the  same 
results  ns  decreasing  frequency  and  vice  versa. 


Starting  from  these  observations,  an  elegant  method  of  reducing  the 


dynstnic  behavior  of  elastotners  (and  other  polymers)  j  covering  a 
range  of  frequencies  and  temperatures,  to  a  single  curve  relating 
appropriate  reduced  properties,  was  developed  by  Ferry,  Fitzgerald, 
Grandine  and  Williams  (ref.  13).  The  method  is  sometimes  referred 
to  as  the  WLF  law,  due  to  a  paper  by  Williams,  Landell  and  Ferry 
(ref.  14),  in  which  the  mechanics  of  its  implementation  are  described 
and  demonstrated.  The  mathematical  statement  of  this  method  is  as 
follows 


(22) 

— 8  . 86  (9  9  ) 

c 

101 .6  +  (0—0  ) 
c 

(23) 

9 

■  <‘=><0.0  f 

(24) 

^  characteristic  temperature  given  by 

^  =  ®t  (50°  ±  5°)  K  (25) 

9^  is  the  transition  temperature  between  the  rubbery 
and  glassy  states  of  the  material. 

G  is  the  observed  modulus  at  temperature  9  and 
frequency  o). 
is  the  reduced  modulus, 

Bq  is  the  coefficient  by  which  test  frequency  is 

multiplied  to  obtain  the  corresponding  reduced 

frequency  at  temperature  0  . 

c 

By  these  relationships,  knowledge  of  a  modulus  as  a  function  of  fre¬ 
quency  at  several  temperatures  may  be  interpreted  as  knowledge  of 
reduced  modulus  over  a  wide  frequency  range  at  the  characteristic 
temperature.  This  knowledge  may  be  further  interpreted  to  give  the 
variation  of  modulus  over  a  wide  frequency  range  at  any  temperature 
between  9^  and  the  maximum  temperature  tested. 

The  method  has  been  demonstrated  to  be  applicable  to  material  prop¬ 
erties  of  storage  and  loss  moduli  or  to  element  properties  of  stiff¬ 
ness  and  damping.  The  additional  information  needed  to  implement 
and  use  the  methods,  beyond  the  elastomer  dynamic  data  is  the  value 

of  9  . 
c 

The  transition  temperature,  0^,  may  be  measured  independently  or  0 
may  be  inferred  from  a  sufficiently  large  volume  of  dynamic  data.  ^ 


where 


and 


i 

ni 


S 

d 


m 


24 


Inf  luence  oi  large  amplitude!  displacements.  Observed  properties 
of  an  (*iastonicr  element  tend  to  be  u  function  of  amplitude  of 
loading,  as  well  as  frequency  and  temperature.  This  is  shown,  for 
ex^niiple,  by  references  15,  16  and  17.  No  specific  and  complete 
explanation  of  this  mechanism  has  been  produced.  It  may  be 
hypothesized,  however,  that  a  major  contribution  is  the  increase 
in  temperature  of  the  element  when  a  significant  dissipation  rate 
is  imposed.  As  will  be  discussed,  further,  in  relation  to  the  test 
data  a  second,  nonthermal,  mechanism  has  also  been  shown  to  cause 
degradation  in  stiffness  at  high  strains. 


INTERPRETATION  OF  BASE  EXCITATION  RESONANT  MASS  TEST  DATA 


The  preceding  review  of  elastomer  material  properties  was  presented 
as  a  framework  against  which  the  methods  of  dynamic  testing  described 

in  this  report  and  the  methods  of  data  reduction  and  interpretation 
may  be  viewed . 


The  BERM  method  has  been  applied  in  the  generation  of  dynamic  prop 
erties  of  three  categories  of  elastomer  element:  shear  elements- 
compression  elements;  and  ring  cartridge  elements.  The  shear  and 
compression  elements  have  been  illustrated  in  Figure  3.  The  ring 
cartridge  element  is  shown  in  Figure  10,  '  ^ 


Table  I  summarizes  the  range  of  conditions  covered  by  this  data  and 
the  periods  when  it  was  generated.  The  1973,  Phase  II,  data  was 
previously  reported  by  Gupta,  Tessarzik  and  Cziglenyi  (ref.  2)  in 
terms  of  stiffness  and  damping  as  a  function  of  frequency.  To 
obtain  a  continuous  mathematical  representation  of  the  dependence 
of  stiffness  and  damping  on  frequency,  a  power  law  relationship  was 
fitted  to  the  test  data  and  the  results  of  this  data  analysis  are 
presented  in  Table  II. 


Prior  to  the  generation  of  the  Phase  III  test  data  presented  in  this 
report,  an  effort  was  undertaken  to  develop  methods  of  exploiting 
the  information  generated  by  the  BERM  test  method,  to  yield  material 
properties  and  geometry  factors  which  can  be  applied  in  engineering 
esign.  The  data  of  Table  II  was  used  to  implement  and  test  these 


Generation  of  Material  Properties 

In  the  previous  section,  it  was  shown  that  stiffness  and  damping 
values  can  be  interpreted  in  terms  of  relationships  between  stress 
and  strain.  It  was  also  shown  that  combined  stress  effects  can  be 
instrumental  in  making  the  effective  moduli  for  the  shear  and  com¬ 
pression  element  different  from  uniaxial  moduli.  Table  III  sum¬ 
marizes  the  static  stiffness  correction  factor  (Ggff)  IQ'  which 
would  be  expected  on  the  basis  of  Eqs .  (12)  and  (13)  for°the  shear 
and  compression  elements  tested  under  Phase  II  and  III. 


In  all  but  two  of  the  elements  tested  the  predicted  effects  of 
combined  stress  are  greater  than  ten  percent.  While  this  data  is 
for  static  behavior  of  the  elements,  and  no  rigorous  basis  exists 
for  applying  it  to  dynamic  behavior,  it  is  assumed  that  the  two 
shear  elements  for  which  the  static  stiffness  correction  factor  is 
within  two  percent  of  unity  are  also  likely  to  show  small  combined 
stress  effects  under  dynamic  loading  conditions.  AIL  the  other 
eioments  tested  are  Likely  to  show  significant  combined  stress 
effects , 
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TABLE  I 

SUMMARY  OF  DYNAMIC  TEST  CONDITIONS 
Material  Tested  -  Polybutacliene 


Dimensions  (eni) 


Phase  Date  Geometry 


II  1973 


III  1975  Shear 


7-4  Compression 


r 

h 

X 

b  d 

0.3175 

0.508 

4.877 

0.3175 

1.27 

4.877 

0.3175 

2.54 

4.877 

0.3175 

1.27 

0.635 

1.27 

1.27 

1.27 

0.3175 

1.27 

4.877 

0.3175 

1.27 

4.877 

0.3175 

1.27 

4.877 

0.3175 

1.27 

4.877 

0.3175 

2.54 

4.877 

0.3175 

2.54 

4.877 

0.3175 

2.54 

4.877 

0.3175 

2.54 

4.877 

0.3175 

2.54 

4.877 

0.3175 

2.54 

4.877 

0.635 

1.27 

0.635 

1.27 

0.635 

1  .27 

0.635 

1  .27 

0.635 

1  .  27 

0.635 

.1.27 

Cartridge 


Nomi na  1 

N 

Temper¬ 

ature 

Dissipation 
W /cm ' 

20 

32  C 

0.022 

8 

4 

32  C 
32  C 

0.022 

0.022 

3 

32  C 

0.022 

10 

30 

32  C 
32  C 

0.022 

0.022 

8 

32  C 

0.062 

8 

50  C 

0.062 

8 

66  C 

0.06  2 

3 

93  C 

0  .  C>  6  2 

4 

32  C 

0.033 

4 

32  C 

0.  176 

4 

32  C 

0.  352 

4 

32  C 

0.704 

4 

4 

32  C 

32  C 

!  .  76 

4.30 

10 

5  C 

0.062 

10 

66  C 

!  0.062 

10 

93  C 

0.062 

10 

32  C 

0.43 

10 

32  C 

1  .  1 

10 

32  C 

1  .905 

2  .858 

0.476 

'4 

32  C 

0.022 

1.905 

2  .858 

0.476 

4 

43  C 

0.022 

1.905 

2.858 

0.476 

4 

66  C 

0,022 

1  .905 

2  .858 

0.476 

4 

80  C 

0.022 

1.905 

2.858 

0.476 

4 

32  C 

0.11 

1 .905 

2  .858 

0.476 

4 

32  C 

0.22 

1.905 

2  .858 

0.476 

4 

32  C 

0.88 

1.905 

2.858 

0.953 

4 

32  C 

0 . 0  3  3 

1.905 

2.858 

0.953 

4 

43  C 

0.03  4 

1.905 

2  .858 

0.953 

4 

66  C 

0.03  3 

1.905 

2  .858 

0.953 

4 

80  C 

n  n./.  / 

TABLE  II 


SUMMARY  OF  ELASTOMER  TEST  DATA  FROM  PM3E  II  TESTS 


MATERIAL  «  Pol ybutadiene 

TEMPERATURE  *  32  deg-C 

CUMPLEX  STIFHNESS  K  •  Kl  ♦  1K2  Ib/in. 

CORRELATIUN  FORM  Y  *  A*OMEGA**B  TH  rad/sec) 


SPECIMEN  ★**♦*****★★*★*  Kt  #**♦**★♦♦  #★♦*&★** ir^**^*  K2  **********a^**  i^****#**#**  AMPLITUDE 


DETAILS 

MEAN 

A 

B 

MEAN 

A 

B 

MEAN 

A 

Ib/in 

lb/ In 

in 

0.50IN-30/C/2.5 

/.54IE*04 

2.224E*04 

I.589E--01 

8.347E'^C3 

5  .873E+04 

-1 .051E“0S 

1 .055E-03 

3.854E-0i  -7.67f^- 

0.50IN-30/C/5.0 

7.809E*04 

2.545E*04 

1.450E-OI 

8.657E+03 

1  .507E^C'! 

-7.2I2E-02 

U042E-03 

3.401 E-OI  -/o52o  - 

0.25IN-IO/C/2.5 

a.042E*04 

2.871E+04 

l.325E»01 

9.2d6E+03 

2.809E+04 

-1.427E“0S 

3.806E-04 

2.040E-01  -8.062 

0,25IN-l0/C/5.0 

8.026E+04 

3.577E+04 

1 .044&--0I 

9.447E+03 

2.753E+04 

-i .3S2E-0I 

3.924E-04 

i,764E-0l  -?.gav;- 

0.I251N- 

3/C/2.B 

7.662E+04 

7.312E+04 

6.368E-03 

1 . J90£*04 

I.378E+04 

»2.004E-02 

2.501 E-04 

2.407E-02  -6.256(7- 

0.)25IN- 

3/C/5.0 

d.392E4'04 

7.007E+04 

3.227E-02 

1 .295E+04 

1 .054E*O4 

-4.859E-02 

2.250E-04 

3.547E-02  -6.854E- 

0.51 N>8 

/S/0. 

1 . 790E*O5 

5.330E*04 

I.548E-0I 

1 . 720E+04 

3.457E*04 

-8.885E-02 

3. 106E-O4 

2.478E-01  -8.473E- 

0.5IM-8 

/S/2.S 

1 . 666E*05 

5.003E^04 

I.530E-01 

l•737E*04 

3.I39E+04 

-7.521 E-02 

3.49!?E-04 

4.985E-0!  -9.232E- 

0.61M-8 

/S/5.0 

1 . 598E*05 

5.36)E*04 

I.39IE-01 

1 .653E+04 

2.678E*04 

-6,l42E-02 

3.462^-04 

2.162E-01  -8.198S- 

I.OlN-4 

/S/0. 

1  .453E+05 

2.569E*04 

2.  I 99E-0I 

1 .624E+04 

2.800E*04 

-6.9I2E-02 

3.497E-.  4 

2*?44E-0)  -a.503E- 

1 .OlN-4 

/5/2.5 

U35IE+05 

2.550E+04 

2.I35E-01 

1 .63IE+04 

2.388E+04 

-4.882E-02 

3.799E-U4 

-8.309E- 

1.0  IN-4 

/S/5.0 

! . 307E*05 

24263E+04 

2.230E-OI 

I .563E+04 

1.923E*04 

-2.647E-02 

3.820E-04 

-a.339E- 

0.2  IN-20 

/S/0. 

).625£*05 

3.708E*04 

I.88IE-0I 

^^89E*04 

5.700E+04 

-i .548E-0! 

3.395E-04 

:'.i6eE-os  -a,220E- 

0.21N-20 

/S/2.5 

1 .549E*05 

3.33dE*04 

1 .936E-0I 

1 .68IE+04 

5.867E+04 

-I.576E-0I 

3.398E-04 

2.O.30E-V  '  8.386E- 

0.2IN-20 

/S/5-0 

! .452E*05 

3.689E+04 

1 .743E-0I 

1  .600E«04 

4.833E*04 

“1  .406E-0  5 

3.683E-04 

2.240E-0i  -8.1562- 

CUMPLEX  STIFFNESS  K  «  Kl  ♦  iK2  New/m 
CORRELATIUN  FORM  Y  *  A*OMEGA**B 


SPECIMEN  ★★★★★★**★★★★★*  Kl  «***★*★★★*****  ******** K2  **************  ***********  AMPLITUDE  ******** 

details  mean  a  B  mean  a  B  MEAN  A  B 

New/m  New/m  m 


0.50 IN-30/ c/2. 5 

1 . 32 I E*07 

3.894b*06 

l.589b-0l 

1 .462E+06 

3.280E*06 

-1.05IE-01 

2.679E-05 

9,790E-03  -7.677E- 

0.50IN-30/C/5.0 

I.368E+07 

4.457E+06 

1.458E-0I 

1 .5I6E+06 

2.640E+06 

-7.212E-02 

2.646E-05 

8.639E-03  -7.529E- 

0.25IN-10/C/2.5 

1 .4O0E*O7 

5.027E+06 

1 . 325E-0 1 

1 .623E+06 

4.920E+06 

-1 .427E-01 

9.668E-06 

5.180E-03  -a.082E- 

0.25IN-I 

O/C/5.0 

1 .406E*07 

6.264E+06 

I.044E-01 

1 .654E*06 

4.022E*O6 

-1 .seaE-ci 

9.968E-06 

4.40OE-O3  -7.a89E- 

0. 125  IN- 

3/C/2.5 

1 -342E407 

1  .281E+07 

6.370E-03 

2.0e4E*06 

2.4I4F*06 

-2.004E-02 

6.354E-06 

6.  M4E-04  -6.216E- 

0.1  25 IN- 

3/C/5.0 

l.557E*07 

1 .227E*07 

3.227E-02 

2.268E+06 

3.247E*06 

-4.059E-O2 

5.7I5E-06 

9.010E-04  -6.8548- 

0.51N-8 

/S/0. 

3. 149E*07 

9.335E*0d 

I.548E-0I 

3.0l2b+06 

6.054E+06 

-a.806E-O2 

8.092E-06 

6.294E-03  -8.4738- 

0.5IN-8 

/S/2.5 

2.9I8E*07 

8. 762E+06 

).530E-0I 

3.042E*06 

5.497E+06 

-7.521E-02 

8.80bE-O6 

I.266E-02  -9,232E- 

0.51N-d 

/S/5.0 

2.799E*07 

9.389b+06 

1.39IE-OJ 

2.095E*Od 

4.689E+06 

-d.l41E-02 

8.794E-06 

5.491 £-03  -8. 190E- 

1 .OlN-4 

/S/0. 

2.545E+07 

4,4  99E+06 

2.I99E-01 

2.844E*06 

4.9038*06 

-6.9i2E-02 

8.08IE-O6 

7.223E-03  -8.503E- 

I.OlN-4 

/s/2.5 

2.366E*07 

4.465E+06 

2.I35E-0I 

2.057E*O6 

4.I83E+06 

-4.883E-02 

9.64bE-06 

6.356E-03  -a.309E- 

1 .0  IN-4 

/S/5.0 

2.289E*07 

3.963E*06 

2.238E-0I 

2.737E+06 

3.369E+06 

-2.647E-02 

9.703E-06 

6.692E-03  -8.339E- 

0.2  IN-20 

/S/0. 

2.846E*07 

6.494E*06 

1 .68IE-0I 

2.959E+06 

9.904E^O6 

-1  .548E-0I 

8.624E-06 

5.506E-03  -8.220E- 

0.2IN-20 

/S/2.5 

2.713E+07 

5.843E+06 

I.936E-01 

2,944E*06 

J .028E+07 

-I.576E-01 

8.631E-06 

6.680E-03  -8.386E- 

0.21N-20 

/S/5.0 

2.542E*07 

6.461 E*06 

I .743b-0l 

2.803E+06 

8.464E*06 

-1 .4068-01 

9.355E-06 

5.691 E-03  -8. I56E- 

N.H.  Frequency  UMFGA  in  the  correlation  function  H  In  red/sec. 


A  sample  consists  of  several  specimens  of  elastomer  arranged  in  a  prescribed 
configuration. 

Specimen  details"  code  given  in  the  table  is  derived  as  follows: 

(1)  C-ompress  i  on  Sample:  (Specimen  thickness)  IN  -  No.  of  Specimens/C/ 

(7o  Preload). 

The  diameter  of  all  compression  specimens  is  1/2 

(2)  Shear  Sample:  (Specimen  length)  IN  -  No.  of  Specimens/S/ 

(7o  Preload) 

The  specimen  thickness  and  width  normal  to 
direction  of  deformation  for  all  shear  specimens 
are  respectively  1/8  in  and  1.92  in. 

29 


TABLE  III.  ELASTOMER  TEST  ELEMENTS  -  EFFECTS 
OF  COMBINED  STRESSES 


Element  Type 

Length 

(cm) 

No. 

Diameter 

(d,cm) 

Strained 

Dimension 

(h,cm) 

Width 

(b,cm) 

^^eff^o 

G 

0 

Shear 

2.54 

4 

0.3175 

4.877 

0.995 

Shear 

1.27 

8 

0.3175 

4.877 

0.980 

Shear 

0.508 

20 

0.3175 

4.877 

0.884 

Compression 

30 

1.27 

1.27 

1.125 

Compression 

10 

1.27 

0.635 

1.5 

Compression 

3 

1.27 

0.3175 

3.0 

This  assumption  leads  to  the  following  approach  to 
the  tost  data. 

interpretation 

1.  From  dynamic  stiffness  and  damping  data  for  the  shear 
elements  of  1,27  cm  and  2,54  cm  length  (parallel  to 
loading  direction)  the  material  properties  of  storage 
and  loss  modulus  should  be  inferred  -  neglecting  what¬ 
ever  small  effects  of  geometry  occur. 

From  these  elements,  values  of  storage  and 
at  a  particular  frequency  are  calculated  as 

loss  modulus 
follows : 

Ki  ((t>)  h 

G'(cd)  =  - 

N  A 

(26) 

K„  (cn)  h 

g"(od)  =  - 

N  A 

(27) 

With  values  of  stiffness  and  damping  available  for  more 
than  one  specimen  at  a  particular  frequency,  the  values  of 
0  (tt) ,  G  (cd)  are  calculated  which  minimize  the  rms  deviation 
between  stiffness  and  damping  predicted  using  the  G'(ai), 

G  (a>)  and  the  measured  stiffness  value  for  the  element. 


to 

to 

to 


2.  From  test  data  for  the  compression  elements  of  different 
shape  should  be  extracted  empirical  coefficients  which 
account  for  the  effects  of  geometry. 

With  material  property  data  available  from  shear  specimen  tests, 
values  of  shape  factor  coefficient  are  determined  which  minimize 
the  deviation  between  predictions  and  measurements  of  compres¬ 
sion  specimen  stiffness  and  damping.  Separate  coefficient  values, 
p  and  p  ,  are  determined  for  stiffness  and  damping  so  that 


Kj^(cd) 

=  G'(co) 

TT  D  N 

4h 

(2:' 

K^Cco) 

=  G'V) 

TT  D^N 

4h 

[i  +  B*  s^] 

(29) 

Implementation  of  Methods  of  Data  Interpretation 

to  Using  the  Phase  II  shear  element  data  for  three  preloads  and  two 

length  values  (six  sets  in  all)  rms  best  fit  values  of  g'  and  g" 
were  determined  over  the  frequency  range  100  to  1000  Hz  and  the 
ig  resultant  set  of  discrete  G  ,  G  frequency  points  were  then  refitted 

to  power  law  expressions.  As  a  result: 


to 


to 


G'(cd)  =  4.136  X  10^ 


(30) 


G  (cd)  =  3.369  X  10^  (31) 

Using  the  Phase  II  compression  element  data  for  two  preloads  and 
three  different  heights,  independent,  best  fit  values  of  p",  p'"  were 
determined  at  discrete  frequency  values  between  100  and  1000  Hz  and 

their  frequency  dependence  was  then  expressed  as  a  power  law  relation¬ 
ship  so  that 


GgffCtu)  =  G'(a))  1^1  +  12.33  od"0-2896 
Ggff(a))  =  g"(cd)  [l  +  1.726  (^)^ " 


(32) 

(33) 


Verification  of  Methods  For  Determining  Material  Properties  And 
Geometry  Factors  . 

Based  upon  the  preceding  results,  individual  shear  elements  of  poly¬ 
butadiene  would  be  expected  to  exhibit  the  following  values  of  stiff¬ 
ness  and  damping  (neglecting  bending  effects): 


to 
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K  (oo)  =  4.136  X  10^  ^  (34) 


K„  (cd)  =  3.369  X  10^ 


(35) 


and  compression  elements  of  polybutadiene  would  be  expected  to  exhibit 
the  following  values  of  stiffness  and  damping: 


K-  (03)  =  4.136  X  10  03 


6  0 . 184  TT  D 


2  - 


4h  L 


1  +  12.33  03 


-0,2896 


4h;  J 


(36) 


V  fn-.\  -  T  -ifLQ  1  -0.0793  ttD 

K„  (03)  =  3  .369  X  10  03  7; — 

^  4h 

To  determine  the  "goodness  of  fit"  given  by  the  above  expressions 
with  the  measured  data  from  which  they  were  derived,  predicted  stiff¬ 
ness  and  damping  of  shear  and  compression  specimens  are  presented  in 
Figures  11  through  22.  Three  different  frequencies,  100,  300  and 
1000  Hz  have  been  selected  and  at  each  of  these,  stiffness  and  damp¬ 
ing  has  been  plotted  against  a  geometrical  parameter  representative 
of  the  difference  among  the  specimens  of  a  particular  class.  For 
the  shear  specimens,  this  parameter  is  the  ratio  of  sheared  length 
to  thickness  (X/xh),  for  which  data  over  a  range  from  0.2  to  1.0  is 
available.  For  the  compression  specimens,  the  independent  variable 
is  the  ratio  of  diameter  to  height  (D/h),  for  which  data  over  a  range 
from  1  to  4  is  available. 

The  consistency  between  the  predictions  and  the  test  data  is  generally 
reasonable,  although  scatter  in  the  test  data  leads  to  individual 
discrepancies  of  up  to  20  percent  of  the  predicted  value.  The  mean 
discrepancy  is  about  ten  percent  of  the  predicted  value. 

The  results  presented  in  Figures  11  through  22  are  taken  as  providing  a 
satisfactory  check  of  the  methods  presented  for  modeling  and  data 
interpretation.  The  same  approach  to  development  of  material  property 
data  from  shear  specimen  tests  will  be  applied  to  the  Phase  III  data 
to  be  presented  subsequently.  This  data  covers  a  much  wider  temper¬ 
ature  range  and  also  reflects  improvements  in  the  test  method  which 
evolved  during  the  Phase  III  tests. 

Alternative  Correlations  For  Shape 

The  work  of  Nashif,  Cannon  and  Jones  (ref.  10)  provides  data  on  the 
stiffness  of  cylinders  in  compression.  Dynamic  stiffness  values  for 
®  range  of  diameter  to  height  ratios  between  0.25  and  4.85  are  given 
at  250  and  500  Hz.  The  material  tested  is  paracryl  (a  Buna  Nitrile 
elastomer  manufactured  by  UniK’oyal).  Just  as  for  the  polybutadiene 


1  +  1.726  CD 


0.02992 


4h/ 


(37) 
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Fig.  20  Phase  II  Correlation  for  Shape  -  Sinsila  Compression 
Element.  Diameter  0.5  In.,  100  Hz  Damping 


or  Shape  -  Single  Govir  ession 
In.,  ;}00  Hz  Damping 


la t ion  for  Shape  -  Single 
ter  0.5  In.,  1000  Hz  Damp 


Jashif,  Cannon  and  Jencb  (R 
idrlcal  Elastomer  Test  Spec 


Correlation  with  Work  of  Nashif,  Cannon  and  Jones  (Ref, 
Dynamic  Stiffness  of  Cylindrical  Elastomer  Test  Specimen 


stiffness  data,  this  set  of  data  has  been  used  to  calculate  shape 
factor  coefficient  values  for  the  coi'relation: 


Stifiritss: 


(38) 


and  the  results  for  250  and  500  Hz  are  P'  =  1.484  and  1.485,  respec 
tively.  The  corresponding  results  for  polybutadiene  are  1.45  and 
1.29.  Comparisons  between  these  correlations  and  the  test  data  of 
reference  11  are  plotted  in  Figures  23  and  24  in  terms  of  the  ratio 
of  actual  stiffness  to  stiffness  predicted  for  uniaxial  stress. 


This  independent  confirmation  of  the  ability  of  a  square  law  correla¬ 
tion  between  stiffness  correction  factor  and  diameter  to  height  ratio 
to  match  test  data  adds  support  to  the  use  of  such  an  empirical 
approach  for  consolidating  the  present  compression  data. 

Use  of  a  Constant  Shape  Factor  Coefficient 

As  discussed  previously  in  this  report,  the  most  commonly  used  value 
of  static  shape  factor  coefficient  O),  for  cylindrical  compression 
elements  is  two.  A  very  straightforward  method  of  predicting  dynamic 
stiffness  and  damping  could  be  based  on  this  coefficient  value,  as 
follows 


(to) 

=  3G'(a>)Sj^[l  +  2(|;)'] 

(39) 

T^2  -  /  X  2-, 

K2  (to) 

=  3G"«.)  [1  +  2  (i)  J 

(40) 

No  analytical  justification  exists  for  this  application  of  static 
coefficients  but,  since  it  has  the  appeal  of  simplicity  and  generality, 
the  likely  errors  in  the  approach  have  been  tested.  Superimposed  on 
Figures  17  through  22  are  the  above  expressions.  At  high  frequencies 
the  predictions  of  these  expressions  for  stiffness  are  50  percent  high 
relative  to  the  more  representative  lines  based  on  frequency— dependent 
coefficients.  At  lower  frequencies  the  discrepancy  for  stiffness  is  less 
and  the  discrepancy  for  damping  is  closer  to  ten  percent  under  all 
conditions . 

Effects  of  Temperature  And  High  Dissipation 

Since  the  Phase  II  data  did  not  include  temperature  as  a  varied 
parameter  it  cannot  be  used  to  evaluate  any  means  of  representing  the 
effects  of  temperature.  However,  in  anticipation  of  the  Phase  III 
data  lor  high  temperature  and  high  dissipation,  a  model  intended  to 
tie  the  efL'oct.s  of  high  dissipation  to  the  change  in  properties  as  a 


result  of  temperature  Increases,  generated  by  energy  dissipation  in  the 
elastomer,  has  been  developed.  This  model  simply  needs  to  have  avail¬ 
able  values  of  G  and  G  ^  as  a  function  of  frequency  and  temperature 

G'(cd,T)  (41) 

G^Co),!)  (42) 

The  hypothesis  of  this  model  is  that,  when  large  amplitude  vibrations 
are  imposed  upon  an  elastomer  element,  they  give  rise  to  dissipation 
of  energy  throughout  the  element.  As  a  result  there  is  a  tendency 
for  the  elastomer  to  heat  up.  The  temperature  distributions  in  the 
elastomer  are  governed  by  the  equations  of  thermal  conduction  and 
appropriate  boundary  conditions.  As  a  result  of  the  distribution  of 
temperature  in  the  elastomer,  and  the  dependence  of  elastomer  moduli 
on  temperature  there  is  a  distribution  of  material  properties  through 
the  elastomer.  Stress  equilibrium  within  the  elastomer,  accounting 
for  this  distribution  of  properties,  governs  the  forces  generated  by 
the  elastomer,  and  its  apparent  stiffness  and  damping. 

The  implementation  of  such  a  model  in  two  or  three  dimensions  could 
be  very  complicated.  However,  to  typify  the  results  given  by  such  a 
model,  and  to  identify  the  major  items  of  information  needed,  while 
minimizing  the  complexity,  a  one-dimensional  development  has  been 
undertaken . 

An  important,  valid,  assumption  is  that  the  time  scale  associated 
with  the  vibrational  phenomena  (less  than  0.01  second)  is  short 
compared  to  the  time  scale  associated  with  changes  in t he  thermal 
condition  of  the  material  (h^  p  C  /k,  which  is  of  the  order  of  40 
seconds  and  higher).  ^ 

One-Dimensional  Model  For  Large  Dissipation  Effects 

This  model  is  directed  at  either  a  shear  or  a  compression  element  in 
which  uniaxial  stress  prevails.  The  temperature  distribution  is  one¬ 
dimensional  and  the  flow  of  heat  is  normal  to  the  stressed  area.  The 
equations  to  be  solved  are: 

1.  The  energy  equation,  in  terms  of  cycle  average  temperature,  T 


2.  The  dynamic  stress-equilibrium  equation  for  uniaxial  stress 


(44) 


where  U  is  the  amplitude  of  displacement  relative  to  the  base 
of  the  elastomer  element. 


li 


3.  The  dependence  of  properties  on  frequency  and  temperature 

C*  -  G*  (a.,T)  (45) 

4.  Boundary  conditions 

T  =  'fg  at  Z  =  0,  h  (46) 

U  =  0  at  Z  =  0  (47) 

5.  For  imposed  displacement  at  the  top  of  the  element 

U  =  Y*  at  Z  =  h  (48) 

6.  For  base  excitation,  the  equation  of  motion  for  the  supported 
mass 

G'^  ^  A  Y  ~  =  Mcd^Y*  (49) 

where  Y  is  the  amplitude  of  displacement  of  top  relative 

to  base, 

■)V  ^ 

Yq  Is  base  excitation  amplitude  (absolute), 

Y  is  a  factor  equal  to  one  for  shear  elements, 
three  for  compression  elements, 

M  is  the  supported  mass, 

CD  is  the  base  excitation  frequency 
A  is  the  stressed  area,  and 
Tg  is  the  temperature  at  base  and  top. 

The  complete  development  of  Equations  (43)  to  (45)  is  presented  in 
Appendix  A. 

_ the  Equations.  As  a  result  of  the  dependence  of  G*  on 

temperature  the  system  of  equations  is  nonlinear.  While  steady-state 
results  are  of  most  Interest,  solution  in  the  time  domain  has  been 
selected  as  more  convenient  to  implement  than  an  iterative  solution. 
The  solution  method  is  a  modification  of  the  second  order  Crank- 
Nicholson  Method  (ref.  17),  in  which  nonlinear  terms  are  evaluated  by 
extrapolation  forward  one-half  time-step. 

The  energy  equation  is  written  to  apply  at  the  time  t  +  At/2. 
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(41 

"k^l 

1  +  1/2  j|- 

[k  ^  1 1 

IdZ  1 

szJ 

Jj.  '  laz 

L  az  J  j 

t  +  At 


-  ^t> 


=  {-  Y 


^  CO  1 

3Z  2  J 


t+At/2 


VAt/2 


left~liand  side  has  a  linear  operator  and  consists  of  known 
quantities  at  time  t  and  unknown  quantities  at  time  t  +  At.  The 
right-hand  side  is  unknown  since  it  is  a  function  of  temperature  at 
time  t  +  At/2,  but  can  be  approximated  by  the  linear  extrapolation 
from  known  quantities: 

‘‘t+A/2  ■  «/2)E^  -  (l/2)Rj.^^  (51) 

The  equation  for  U  is  solved  at  time  t  +  At,  once  T  and  the 
associated  values  of  G*  are  known.  t+At 

The  time  integration  solution  is  continued  until  results  for  succes¬ 
sive  time  steps  Indicate  that  thermal  equilibrium  has  been  reached. 

The  equilibrium  distribution  of  temperature,  material  properties, 
and  deflection  will  determine  the  behavior  which  will  be  observed 
once  the  initial  thermal  transients  have  died  out. 

An  interesting  feature  of  the  thermal  transient  behavior  occurs  when 
the  base  excitation  boundary  conditions  apply  and  the  frequency  of 
excitation  and  supported  mass  are  such  that  the  element  causes  sub¬ 
resonant  vibrations  when  the  element  is  cold  but  vibrates  above 
resonance  when  the  element  gets  hot  due  to  self-heating.  In  this 
case  there  is  a  wide  swing  in  phase  angle  during  the  thermal  tran¬ 
sient  and  an  inability  to  hold  subresonant  conditions  for  levels  of 
dissipation  higher  than  a  certain  limit.  Figure  25  illustrates  the 
behavior  to  be  expected  for  a  range  of  different  input  amplitudes 
and  this  phenonmenon  is  discussed  further  in  the  Results  and  Discussion 
Sections  of  this  report.  Observed  stiffness  and  damping  of  the  ele¬ 
ment  at  any  time  during  the  time  integration,  including  equilibrium, 
are  calculated  as  * 


Kj  -  ^  Im  {g*  ||}^ 
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High  Dissipation  Predictions  Influence  of  Input 
Amplitude.  1.816  Kg  Mass,  380  Hz,  32  C  Ambient, 
Compression  Specimen,  10  Elements,  1,27  cm 
Diameter,  0.635  cm  High 
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It  is  important  to  note  that  a  number  of  material  properties  not  in¬ 
cluded  in  the  dynamic  tests  are  needed  to  implement  this  analysis;  in 
particular,  density,  specific  heat  and  conductivity. 

In  addition,  if  this  analysis  is  to  be  applied  with  meaning  to  com¬ 
pression  elements  in  which  the  stiffness  correction  factor  is  of  order 
1,5  and  above,  some  account  must  be  taken  of  the  stiffness  correction 
factor.  -In  the  present  case,  the  significant  assumption  has  been 
made  that  the  stiffness  correction  factors  follow  the  functional 
dependence  on  frequency  and  geometry  determined  empirically  for  cylin- 
ricai  elements  and  that  the  increase  in  effective  modulus  should  take 
place  on  a  local  basis,  and  be  reflected  in  the  equilibrium  equation: 


PREDICTION  OF  CARTRIDGE  DYNAMIC  PROPERTIES 
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Oni‘  ol  the  objectives  of  the  work  reported  is  to  develop,  and 
evaluate,  relatively  simple  predictive  methods  for  practical  elastomer 
configurations,  such  as  radially  loaded  ring  cartridges.  In  this 
Section  a  number  of  alternative  methods  of  prediction  of  radial  stiff¬ 
ness  are  developed  and  numerical  values  are  obtained  for  the  particular 
cartridge  configurations  selected  for  testing. 

Figure  10  illustrated  the  geometry  of  the  elastomer  ring  cartridge. 

The  elastomer  material  separates  the  inner  diameter  of  one  cy^'inder 
from  the  outer  diameter  of  a  second,  smaller  cylinder.  The  lo^iling 
of  interest  is  radial;  that  is,  it  tends  to  move  one  cylinder  radi.'ly 
eccentric  to  the  other. 

Three  predictive  methods  have  been  tested: 

•  Beam-Column  Method 

•  Method  of  Gbbel  (ref,  19) 

•  Predictions  Based  Upon  Plane  Stress  Analysis 

Each  of  these  is  a  'static*  prediction  method;  that  is,  it  seeks  a 
ratio  of  element  stiffness  to  storage  modulus  which  is  a  function  of 
the  geometry  only  and  applies  that  ratio  in  the  determination  of 
dynamic  stiffness  and  of  damping.  The  approach  is  not  rigorously 
justifiable,  but  may  be  given  credibility  by  empirical  verification. 

Beam-Column  Method 

In  this  development  the  cartridge  is  assumed  to  be  made  up  of  a 
series  of  noninteracting,  elemental  column/beams ,  whose  variation  of 
thickness  due  to  radius  effects  is  small.  Figure  26  illustrates  the 
cartridge  and  the  Beam-Column  element,  whose  mean  loaded  area  is 

da  =  r  i-d9  (55) 

where  i  is  the  length  of  the  cartridge 

F  is  the  mean  radius  =  (r^  +  r^)l2. 

Consider  the  rigid  inner  member  to  be  displaced  radially  by  a  small 

radial  distance  Y,  relative  to  the  rigid  outer  member.  The  compression 

deflection  of  the  element  is  Y  cos  9,  ii  9  is  the  angle  the  radial  line 

to  the  element  makes  with  the  direction  of  displacement.  The  shear 

deflection  of  the  element  is  Y  sin  9.  So,  considering  the  element  to 

act  independently  of  adjacent  elements,  the  elemental  compression  and 

shear  forces,  df  ,  df  are 
c  s 

df  =  Y  cos  9  E  — r  (56) 

c  (r.  -  r  J 
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df  =  Y  8ln  e  G  - r 

s  (r^  -  rp 


and  the  total  force  exerted  at  the  mean  radius  is 


X*  2  o  * 

77  TTTTT  E  cos  6  dO  +  G  sin*  9  d9 
o  '•’^2  ^1'' 


:=  -l-LiJI  I  E  + 

(r^  -  rp  ^  ®J  -58) 

SO  that,  if  we  assume  E  =  3G  (Poisson's  ratio  v  =  0.5  implied),  th.--. 

E  (r  +  r.) 

K  =  Y  77 — - \  BEAM-COLUMN  RESULT  (59) 

A  slightly  different  result  is  obtained  if  the  tapering  of  beams  and 
columns  due  to  radius  effects  is  accounted  for. 

K  =  GZ  ACCOUNTING  FOR  ... 

r  i  (r^/r,)  RADIUS  EFFECTS 
n  4  1  "  . ■  "  ' 

But,  as  Table  IV  shows,  the  reduction  in  stiffness  due  to  radius  effect 
is  less  than  four  percent  for  radius  ratios  up  to  two. 


TABLE  IV.  REDUCTION  IN  STIFFNESS  DUE  TO  RADIUS  EFFECTS 
(Relative  to  Beam-Column  Result) 


v^/v^ 


Reduction 


0.08' 
1.35' 
3.827. 


The  radius  ratio  selected  for  testing  is  1.5, for  which  radius  effects 
are  only  1.35  percent,  and  have  been  neglected. 

Method  of  Gobel 

The  method  of  Gbbel  initially  follows  a  very  similar  development  to  the 
Beam-Column  method  accounting  for  radius  effects,  with  the  intermediate 
result : 
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However,  GCSbel  then  claims  that,  for  /  _  \  =  1,  the  effective  E-modulus 

is  given  by  '■^2  ^1' 


E 


eff 


6.5  G 


(62) 


and,  in  addition  that,  for  other  length  to  radius  ratios,  a  correction 
factor,  fj^,  should  be  applied,  so  that 


K 

r 


7.5  ni  Gf. 


(63) 


Gobel  presents  fj^  as  a  graphical  function  of  -  r^).  Payne  and 

Scott  (ref.  6)  in  their  book  ENGINEERING  DESIGN  WITH  RUBBER  present  a 
cubic  fit  to  Gobel 's  function: 

f =  1  +  0.0097  -  r^)^  (64) 

which  agrees  closely  with  Gobel's  graphical  presentation.  The  basis 
for  f]^,  or  for  the  statement  E  =  6.5  G  is  not  made  clear  in  Gobel's 
book.  These  are  not  presented  as  empirical  results,  since  no  test 
data  is  given  for  this  geometry.  Gobel's  method  implies  that  Eff  = 

6.5  G  is  an  asymptotic  condition  for  very  short  length  cartridges. 

This  implication  is  surprising. 

One  set  of  static  test  data  for  the  cartridge  configuration  (Adkins 
and  Gent,  ref.  7)  has  been  identified  and  this  data  is  compared  with 
the  prediction  of  the  Beam-Column  method  and  the  method  of  Gobel  in 
Figure  27.  It  shows  good  agreement  at  higher  values  of  X/(r„  —  ri  ) 
between  Gobel  and  test  data  but  at  lower  i/iv  -  r^^),  the  Beam-Column 
method  appears  to  come  closer,  due  to  a  much  sharper  drop  in  stiffness 
of  the  test  element.  This  comparison  suggests  that  Egff  =  6.5  G  is 
not  a  valid  asymptotic  condition  for  short  cartridges. 

Predictions  Based  Upon  Plane  Stress  Analysis 

The  value  of  using  a  plane  stress  elastic  analysis  was  also  investigated. 
In  this  case,  only  a  quarter  of  the  cartridge  was  modeled  and  appropri¬ 
ate  symmetry  and  antisymmetry  boundary  conditions  were  imposed.  The 
quarter  cartridge  was  embedded  in  a  rectangular  finite  difference  mesh 
and  a  finite  displacement  of  the  inner  member  relative  to  the  outer 
member  was  imposed.  An  existing  computer  program  (MTI  Program  No.  444) 
was  used  to  solve  for  the  internal  displacement  field  within  the 
cartridge.  Since  it  is  not  possible  to  employ  a  Poisson's  ratio  of 
0.5  in  this  computer  program,  a  value  of  0.4983  was  used,  together 
with  a  shear  modulus  value  of  1.795  x  10^  N/m^.  This  combination  of 
Poisson' s^ratio  and  shear  modulus  corresponds  to  a  bulk  modulus  of 
5.28  X  10  N/m  ,  which  is  within  the  range  of  bulk  moduli  determined 
by  Rightmire  (ref.  18). 
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RED 


RATIO  OF  LEN6H  TO  RADIAL  EXTENT 

-^  /  (r2-r|) 

I'i;.;.  27  Comparison  of  Prediction  Metliods  for  Cylindrical  Cartridge 
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The  stiffness  is  calculated  by  determining  the  component  of  shear  and 
direct  stress  in  the  direction  opposing  the  displacement  x  and  inte¬ 
grating  over  the  active  cartridge  surface  to  obtain  the  net  force. 

The  stiffness  is,  then,  the  ratio  of  net  force  to  displacement. 

Numerical  Predictions  of  Cartridge  Stiffness  and  Damping 

Each  of  the  predictive  methods  described  for  the  ring  cartridge  has 
been  applied  to  the  geometry  selected  for  testing.  The  results  are 
presented  in  Table  V  in  nondiraensional  form;  that  is,  the  ratio  of 
stiffness  to  the  product  of  shear  modulus  and  length  is  presented. 


TABLE  V.  NONDIMENSIONAL  STIFFNESS  FOR  CARTRIDGES 
SELECTED  FOR  TESTING  (Inner  and  Outer 
Diameters;  1.905  and  2.858  cm) 


Method 

K/GX  For  i  =  0.476  cm 

K/Gi»  For  SL  =  0.953  cm 

Gobel 

58.7 

62.6 

Plane  Stress 

35.1 

35.1 

Beam-Column 

31.4 

31.4 

It  is  seen  that  only  the  Gobel  method  reveals  any  dependence  of  this 
nondimens ional  stiffness  on  length  of  the  cartridge,  and  that  the  level 
of  the  Gobel  predictions  is  substantially  higher  than  that  of  the  Plane 
Stress  and  Beam-Column  predictions , 


The  assumption  necessary  to  apply  the  Gobel  and  Plane  Stress  predictions 
in  the  determination  of  dynamic  stiffness  and  damping  is  that: 


K 

G  SL 


(65) 


where  +  iK2  is  the  complex  impedance  of  the  cartridge  and  G*  = 

G  +  iG  is  the  complex  modulus.  That  is,  the  relative  distribution  of 
stresses  is  independent  of  time  under  sinusoidal  loading  and  displace¬ 
ment  . 


The  simple  Beam-Column  method,  without  accounting  for  radius  effects, 
since  it  assumes  no  interaction  between  elements,  may  be  developed 
directly  using  either  G  or  G’  .  In  other  words,  beyond  the  significant 
assumptions  in  its  development  already  presented,  no  further  assumptions 
are  necessary  in  applying  it  to  dynamic  stiffness  and  damping. 
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Using  these  assumptions  where  appropriate,  the  predictions  of  Table  V 
may  be  applied  to  the  cartridges  tested  by  the  BERM  method*  The  non- 
dimensional  stiffness  from  Table  V  is  simply  multiplied  by  the  length 
of  the  cartridge  and  the  storage  and  loss  modulus  values  at  the  fre¬ 
quency  of  interest. 

The  results  of  this  prediction  approach  are  compared  with  test  results 
later  in  this  report. 


EXPERIMENTAL  DESIGN  AND  DESCRIPTION  OF  TESTS 

A  forced  vibration,  resonant  mass  type  of  apparatus  was  used  to  test 
the  dynamic  properties  of  the  elastomer  samples.  The  apparatus  was 
designed  as  a  base  excitation,  electromagnetic  shaker  driven  mass- 
spring  system  which  could  be  brought  to  resonance  for  a  range  of  dif¬ 
ferently  sized  vibrating  masses  on  top  of  the  elastomer  ’spring*. 
Acquisition  of  data  was,  however,  not  limited  to  the  resonance  con¬ 
dition  of  tlie  system.  In  fact,  data  obtained  at  the  resonant  frequency 
of  each  mass— spring  combination  was  just  one  of  several  data  points 
acquired  at  several  vibration  frequencies  around  resonance,  where 
significant  amplitude  ratios  between  base  and  mass  vibration  existed. 

The  following  sections  present  a  description  of  the  test  rig,  the 
design  of  the  elastomer  test  samples,  the  instrumentation  as  it  was 
used  for  the  recording  of  test  data,  and  the  range  of  test  parameters 
investigated. 

Description  of  Elastomer  Test  Rig 

The  basic  features  of  this  test  rig  have  largely  been  developed  in  the 
course  of  work  conducted  under  previous  phases  of  NASA  Contract  NAS3- 
15334.  Reference  2  provides  a  complete  description  of  test  rig  capabil¬ 
ities  and  an  analysis  of  the  dynamic  characteristics  of  the  test  rig 
components.  Figure  28  is  a  schematic  showing  all  the  components. 

For  the  tests  reported  herein,  not  all  of  the  available  features  of  the 
test  rig  were  used.  In  particular,  for  test  frequencies  above  100  Hz, 
neither  the  mass  support  air  spring  nor  the  guide  spokes  (guide  bearings) 
were  used.  The  simplified  test  rig  consisted  of  a  shake  table  mounted 
base  plate  to  which  the  elastomer  samples  were  bonded,  and  a  second 
plate  bonded  to  the  free  side  of  the  elastomer  samples.  A  piston  with 
Bellofram  seals  and  various  amounts  of  deadweight  were  attached  to  the 
second  plate  which  vibrated  freely  on  the  elastomer  spring  when  the  base 
plate  was  subjected  to  sinusoidal  vibrations  from  the  shake  table. 

The  piston  cavity  was  externally  pressurized  with  air  and  served  to 
preload  the  elastomer  sample  in  a  static  sense.  An  overall  view  of 
the  simplified  test  rig  preload  cylinder  is  given  in  Figure  29,  where 
a  small  mass  is  attached  to  the  free  end  of  the  preload  piston  pro¬ 
truding  from  the  cylinder.  Preload  was  provided  on  all  tests  conducted 
with  elastomer  compression  samples  and  on  a  few  tests  of  the  cartridge 
configuration.  (The  presence  or  absence  of  preload  is  noted  on  all 
test  results.)  All  tests  with  elastomer  samples  stressed  in  shear  and 
the  majority  of  the  cartridge  tests  were  made  without  static  preload 
on  the  sample,  and  for  these  tests,  the  rig  was  further  simplified  by 
removing  the  piston  and  seals.  Without  piston  and  seal,  all  test  rig 
components  and  features  with  interactive  or  additive  dynamic  contribu¬ 
tions  to  the  test  system  had  been  eliminated.  What  remained  was  a 
dynamically  "clean”  test  rig  which,  of  its  own,  did  not  contribute  any 
identifiable  dynamic  "noise"  or  distortion  of  test  results. 


MIS.' 
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BELLOFRAM  SEALS 


MASS  SUPPORT 
AIR  SPRING 


Cylinder  and  Small  Maas 


The  vibration  input  to  the  test  rig  in  all  tests  was  obtained  from 
a  commercially-available  electromagnetic  ahaker  system^  capable  of 
delivering  66,700  N  (15,000  lb)  force  in  the  sinusoidal  mode  of 
vibration . 

The  ambient  temperature  within  the  elastomer  test  specimen  cavity  was 
controlled  through  the  combined  effect  of  electric  heating  cables  and 
water-cooled  tubes  installed  in  the  base  plate  below  the  test  specimens 
and  heated  or  cooled  air  blown  directly  to  the  test  specimen  cavity. 

For  higher  than  ambient  temperatures,  the  electric  heaters  were  auto¬ 
matically  controlled  by  a  temperature  controller,  which  obtained  the 
necessary  feedback  from  the  thermocouple  embedded  in  the  meta'  of  the 
base  plate.  For  ambient  temperatures  down  to  2  C,  the  air  sup^.  .y  tc 
the  test  specimen  cavity  was  passed  through  cooling  coils  embedded 
dry  ice. 

Details  of  the  temperature  control  system  and  of  the  temperature 
recording  system  are  shown  in  Figure  30, 

Elastomer  Test  Samples  -  Material  Selection 

The  initial  selection  of  polybutadiene  as  test  material  had  already 
taken  place  during  an  earlier  phase  of  this  program,  reported  in 
reference  2.  Restated  briefly,  the  choice  of  polybutadiene,  which  is 
classed  as  Broad  Temperature  Range  (BTR)  elastomer,  was  primarily  based 
upon  the  desire  to  test  a  material  that  was  not  overly  sensitive  to 
small  variations  in  ambient  temperatures  nor  to  the  temperature  varia¬ 
tion  induced  in  the  sample  through  vibration  testing. 

The  choice  of  golybutadiene  among  BTR  elastomers  was  based  upon  the 
manufacturer’s’  recommendation  founded  on  manufacturing  considera¬ 
tions,  such  as  the  achievement  of  consistent  material  properties  from 
batch-to-batch .  This  material,  which  carries  the  manufacturer’s  desig¬ 
nation  NEX156G,  has  a  nominal  hardness  of  70  durometers  (or  Shore  A 
Hardness)  and  has  the  highest  hardness  and  damping  of  those  available 
from  the  manufacturer . 

Elastomer  Test  Samples  -  Configurations 
Three  sample  configurations  were  tested: 

•  Shear  Type  Samples 

•  Compression  Type  Samples 

•  Cartridge  Type  Samples 


Ling  Electronics,  Model  335A  Shaker  with  PP-35/70  VC  Power  Amplifier 
SCO-100  Servo  Control  Center. 


Nichols 


Engineering, 


Incorporated,  Shelton,  Connecticut . 
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Temperature  Control  System  For  Elastomer  Tests 
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S_hear  Test  Samples.  Two  shear  test  samples  were  used  for  testing  in 
this  phase  of  the  contract.  One  shear  test  sample  was  made  up  from 
eight  individual  elastomer  strips,  each  measuring  12.7  mm  in  height, 

48.8  mm  in  length  and  3.2  mm  in  thickness  (0,5  x  1.42  x  0.125  inch) 

(Fig.  31).  A  second  shear  sample  consisted  of  four  individual  elastomer 
strips,  each  having  twice  the  height  (25.4  mm)  of  the  eight-strip 
sample  (Fig.  32). 

One  individual  strip  was  glued,  on  its  large  face  to  each  of  the  sides 
of  a  square  aluminum  block,  with  the  other  large  face  of  each  strip 
glued  to  an  additional  aluminum  holding  block.  A  cyanoacrylate  gine 
was  used  (Eastman  910  MHP) .  In  the  eight-strip  sample,  two  such  assem¬ 
blies,  each  consisting  of  four  strips,  were  placed  on  top  of  each 
other  and  tested  together.  In  the  assembly  of  the  test  fixture,  the 
four  holding  blocks  were  bolted  to  the  shake  table,  while  the  elastomer 
supported  center  block  was  attached  to,  and  therefore  became  part  of, 
the  resonant  mass. 

Compression  Test  Sample.  The  compression  test  sample  consisted  of 
ten  individual  elastomer  cylinders,  each  measuring  12.7  mm  in  diameter 
and  6.4  mm  in  height  (0.5  by  0.25  inch).  These  cylinders  were  arranged 
in  a  circular  pattern  and  glued  with  their  faces  to  two  aluminum  plates 
(see  Fig  33).  The  bottom  plate  was  anchored  to  the  shake  table  and  the 
top  plate  connected  to  the  resonant  mass. 

Like  the  shear  test  sample,  the  compression  test  sample  had  originally 
been  made  for  the  investigations  reported  in  reference  2.  Details  of 
manufacture  and  test  results  for  two  other  compression  test  samples 
made  up  of  three  cylinders  3.2  mm  in  height  and  of  thirty  cylinders 
each  12.7  mm  in  height  may  also  be  found  in  reference  2. 

Cartridge  Test  Samples.  Individual  cartridge  test  specimens  were  made 
up  of  an  annular  elastomer  body  with  rectangular  cross-section  held 
between  the  cylindrical  bore  of  a  rectangular  shell  structure  and  a 
round  shaft  of  smaller  diameter. 

The  elastomer  ring  was  molded  to  the  shaft  and  bonded  to  the  surround¬ 
ing  shell,  which  was  of  split  design  to  facilitate  initial  assembly  and 
later,  static  preloading  prior  to  testing. 

Figure  34  shows  one  cartridge  specimen  before  and  after  assembly  of  the 
outer  shell  around  the  bonded  elastomer  ring  on  the  shaft. 

Cartridge  dimensions  were  chosen  so  that  four  individual  cartridges  in 
parallel  would  have  a  combined  spring  stiffness  that,  in  combination 
with  the  minimum  mass  arrangement  in  the  test  rig  (top  plate  only), 
would  result  in  a  system  resonant  frequency  of  at  least  700  Hz.  Fre¬ 
quencies  to  1000  Hz  could  then  be  achieved  by  testing  to  about  1.4  times 
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MTI-16652 


Test  Assembly  of  Eight  Elastomer  Shear  Specimens, 
Each  1.27  cm  (0.5  In.)  High 
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MTI-16662 


Fig.  34  Elastomer  Cartridge  Specimen  Before  Bonding  Outer 
Shell  (Left)  and  with  Outer  Shell  in  Place  (  ight) 


the  highest  resonant  frequency  and,  through  mass  addition 
frequencies  lower  than  200  Hz  could  be  achieved. 


resonant 


""k  cartridge  specimens  was  designed  with  identical  metal 
dimensions  but  with  only  half  the  elastomer  width  of  the  first  set. 

or  this  set  the  upper  limit  of  the  frequency  range  was  reduced  In 
proportion  to  the  square  root  of  the  stiffness  change,  or  by  a  factor 
of  1.4. 


The  elastomer  dimensions  (in  mm)  of  tlie  individual 
two  sets  were  as  follows: 


specimens  for  the 


Width 

9.5 

4.7 


Outer 

Diameter 

28.6 

28.6 


Inner 

Diameter 

19.1 

19.1 


elastomer  specimens  is  shown 

in  Figure  35.  The  manufacturing  process  of  the  cartridge  elements 
evolved  from  discussions  with  the  manufacturer.  Since  a  split  shell 
was  desired  to  allow  a  positive  radial  preload  to  be  appliL  for 
testing,  molding  of  the  elastomer  to  both  shaft  and  outer  shell  was 
not  considered  feasible  because  tensile  stresses  setup  in  the  elastomer 
ue  to  material  shrinkage  would  not  be  circumferentially  symmetrical. 

"  selected  manufacturing  process,  the  elastomer  was  first  molded 
on  the  shaft  and  then  ground  on  the  outer  diameter  to  the  dimension  of 
the  circular  hole  in  the  assembled  shell.  With  a  small  interference 
o  erance  tf)e  outer  sliell  was  then  clamped  around  the  elastomer  and 

.cilitaud  by  a  special  coating  applied  to  the  bore  of  the  shell  prior 
L  assem  y.  (  le  coating  in  the  shell  is  clearly  visible  in  Fig.  34.) 

Prior  to  assembly  into  a  complete  test  sample,  all  cartridge  elements 
received  a  uniform  preload  of  5.3  percent  of  their  uncompressed  raTi^l 

inv^T*  accomplished  by  replacing  the  shims  measur- 

ing  1.625  mm  between  the  halves  of  the  outer  shell  with  reduced 
hickness  shims  only  1.120  mm  thick.  The  assembly  of  the  Individual 
elastomer  cartridges  into  the  complete  test  sample  is  depicted  in 
Figures  36  through  38.  An  assembly  of  four  cartridges  was  the 
largest  number  which  could  effectively  be  packaged  within  the  exist¬ 
ing  size  constraints  of  the  test  rig.  In  Figure  36  the  four 
elastomer  cartridges  are  shown  next  to  the  top  plate  to  which  they 
will  be  clamped  with  two  clamps  each.  In  the  following  picture 
(Figure  37)  the  cartridges  have  been  placed  into  their  locations 

"fr  ?®’,  ''"■  cartridge  element 

Shafts  and  locked  to  the  plate  with  bolts  from  beneath.  At  this  stage 

le  elevations  of  the  Individual  rnrtri.Ige  element  shells  were  measured. 
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and  shims  were  prepared  for  placement  between  the  flat  surface  of  the 
cartridge  shell  and  the  base  plate  (shown  standing  behind  the  assembly 
in  Fig.  38).  With  shimming  in  place,  differences  in  elevations  between 
individual  cartridge  shells  were  held  to  less  than  0.025  mm.  The 
assembly  of  the  second  set  of  cartridge  elements  (narrow  elastomer 
elements)  received  a  very  light  grinding  cut  to  remove  some  lack  of 
parallelism  between  the  flat  shell  surfaces  and  the  bottom  plate  on 
two  of  the  elements. 

Instrumentation 

The  uncomplicated  measurement  requirements  of  the  base  excitation 
resonant  mass  testing  method  have  been  expounded  before  (see  referej.,  j 
1).  Restated  briefly,  the  measurement  requirements  for  the  experimi."i“ 
tal  investigation  of  the  elastomer  dynamic  properties  are  as  follows: 

1.  Displacement  measurement  of  elastomer  support  plate 
attached  to  the  vibration  table,  relative  to  ground; 

2.  Displacement  measurement  of  elastomer  support  plate 
attached  to  the  resonant  mass,  relative  to  ground; 

3.  Phase  angle  measurement  between  displacement  measure¬ 
ments  (1)  and  (2)  above; 

4.  Displacement  measurement  between  two  elastomer  support 
plates  (relative  amplitude  across  the  elastomer); 

5.  Vibration  frequency; 

6.  Temperature  of  elastomer  support  plates  and  ambient 
air  temperature  in  test  specimen  cavity; 

7.  Vibrational  power  dissipation  across  the  elastomer 
test  sample. 

For  convenience,  acceleration  was  measured  in  (1)  and  (2)  instead  of 
displacement,  but  since  the  motions  were  sinusoidal  the  displacements 
were  directly  derivable  from  the  accelerations.  The  very  low  level  of 
motion  of  the  shaker  body  (from  which  the  shaker  table  is  isolated  by 
air  springs  with  a  very  low  resonant  frequency)  was  also  measured  by 
accelerometers . 

The  input  acceleration  (1)  was  measured  by  an  accelerometer  located  on 
the  support  base  for  the  elastomer  test  samples.  Output  acceleration 
nn  the  resonant  mass  (2)  was  determined  from  three  accelerometers 
mounted  directly  to  the  top  plate  attached  to  the  elastomer  specimens. 

The  use  of  three  accelerometers,  spaced  90  degrees  apart,  with  their 
output  signals  displayed  on  one  oscilloscope  screen,  permitted  immediate 
detection  of  nonaxial  motions  of  the  resonant  mass.  When  nonaxial 
motions  of  the  resonant  mass  occur,  they  manifest  themselves  either  as 
amplitude  or  phase  angle  variations  among  the  three  signals,  or  as  a 
combination  of  both  (testing  experience  revealed  that  noticeable  nonaxial 
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motions  occurred  only  in  a  few  instances  at  certain  test  frequencies. 

Test  data  was  generally  discarded  when  nonaxial  motions  were  observed.) 

The  displacement  measurement  (4)  between  the  two  elastomer  support 
plates  (relative  amplitude  across  the  elastomer)  was  accomplished  with 
a  noncontacting  capacitance-type  sensor. 

Chromel-alumel  type  thermocouples  were  used  to  measure  the  temperature 
in  the  metal  directly  adjacent  to  one  or  two  elastomer  specimens  in 
each  assembly  at  a  distance  of  approximately  1,6  mm  from  the  elastomer 
surfaces.  In  the  shear  sample  the  thermocouples  were  embedded  at  the 
midface  location.  In  the  compression  test  samples  the  thermocouples 
were  directly  above  and  below  the  center  point  of  the  flat  surfaces  of 
the  cylindrical  specimens.  In  the  cartridge  specimens  one  thermocouple 
was  located  in  the  axis  of  the  shaft  at  the  midpoint  of  the  elastomer 
and  tv;o  additional  thermocouples  were  placed  diametrically  opposite  in 
the  metal  of  the  outer  shell,  approximately  1.6  mm  from  the  outer 
diameter  of  the  elastomer  ring  and  also  at  the  axial  midpoint  location 
of  the  elastomer  ring. 

Displacement  and  acceleration  data,  together  with  temperatures  from 
individual  thermocouples  were  displayed  visually  in  digital  form  for 
monitoring  purposes.  As  a  permanent  record,  this  data  was  also 
recorded  on  a  21-column  digital  printer. 

Amplitude  signals,  from  the  four  accelerometers  and  one  displacement 
sensor,  were  sequentially  switched  by  an  analog  scanner  into  a  two- 
channel  tracking  filter,  which  provided  a  visual  readout  of  vibration 
frequency  and  of  two  filtered  amplitude  signals.  The  base  amplitude 
signal  (acceleration)  was  fed  at  all  times  into  one  of  the  two  channels 
of  the  tracking  filter  and  from  there  into  a  phase  meter  where  it  served 
as  a  reference  signal  for  the  measurement  of  the  phase  angle  of  the 
three  output  acceleration  signals  and  the  displacement  signal.  The 
d-c  values  proportional  to  phase  angle  and  amplitude  from  the  phase 
meter  and  the  tracking  filter  were  than  converted  into  digital  form  in 
digital  voltmeters  and  printed.  Each  printed  line  contained  phase  angle 
and  amplitude  from  one  sensor,  together  with  the  vibration  test  fre¬ 
quency  and  the  temperature  indication  from  one  thermocouple.  (See 
Figure  39  for  a  schematic  of  the  data  acquisition  system.) 

To  monitor  the  power  dissipation  in  the  test  sample  during  testing,  a 
previously  developed  power  meter  was  used.  This  meter  utilizes  avail¬ 
able  acceleration  and  displacement  signal  inputs.  The  three  (unflltered) 
accelerometer  signals  from  the  top  plate  of  the  specimen  holder  (resonant 
mass  signals)  are  averaged,  to  give  a  signal  proportional  to  the  net 
force  on  the  resonant  mass.  This  signal  is  therefore  almost  equal  to 
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that  portion  of  the  force  acting  through  the  elastomer  sample  which 
is  due  to  dynamic  displacement  (differing  only  by  dynamic  displace- 
ment  forces  in  the  air  cylinder  in  those  cases  where  static  preload 
is  applied  to  the  test  sample).  The  displacement  signal,  which 
indicates  relative  motion  across  the  sample,  is  differentiated  to 
give  a  signal  proportional  to  rate  of  deformation  (relative  velocity) 
The  force  and  velocity  signals  are  multiplied  together  to  give  the 
instantaneous  power.  The  instantaneous  power  consists  of  a  d-c 
(average)  component  and  a  sinusoidally  varying  component  and,  by 
ilterlng  out  the  sinusoidally  varying  component,  the  average  power 
dissipation  is  obtained. 


Testing  Procedures 

One  of  the  major  features  of  the  current  test  program  was  the  performance 
of  the  scheduled  tests  at  near-resonance  conditions.  Since  phase  angle 
between  base  excitation  and  resonant  mass  response  is  an  accurate  in¬ 
dicator  of  the  amount  of  damping  in  the  region  of  resonance,  measure¬ 
ments  should  preferably  be  made  in  the  phase  angle  range  between  approx¬ 
imately  15  and  165  degrees  (as  discussed  in  reference  2) .  This  requires 
that  the  test  frequency  be  approximately  0.7  to  1,5  times  the  critical 
frequency  of  the  elastomer-resonant  mass  system.  Coverage  of  the  full 
test  frequency  range  of  100  to  1000  Hz  with  sufficient  test  points, 
therefore,  necessitates  changes  in  the  size  of  the  resonant  mass. 
Typically,  between  five  and  seven  different  mass  arrangements  ranging 
from  0.5  to  15  kg  were  found  sufficient  to  cover  the  desired  test  range. 
The  smallest  mass  consisted  of  the  elastomer  test  sample  top  plate  only, 
with  no  additional  weight  attached.  Increasingly  larger  mass  combina¬ 
tions  were  obtained  by  attaching  weights  of  up  to  15  kg  directly  to  the 
top  plate.  In  the  case  of  the  compression  sample,  which  was  tested 
with  preload  applied  through  the  air  piston  cylinder  arrangement, 
weights  up  to  10  kg  were  directly  attached  to  the  preload  piston. 

Once  any  one  of  the  elastomer  test  samples  had  been  mounted  on  the  shake 
table  and  the  necessary  instrumentation  installed  and  connected,  a 
typical  test  sequence  proceeded  as  follows: 

1.  A  resonant  mass  was  selected  and  installed. 

2.  In  the  case  of  the  compression  tests,  the  elastomer 
sample  was  statically  preloaded. 


The  elastomer  teat  sample  cavKy  was  enclosecl  and  the 
temperature  control  system  given  time  to  adjust  ambient 
temperature  to  the  desired  value. 

With  low  vibration  levels  applied  to  the  base  of  the 
elastomer  holding  fixture,  frequency  scans  were  con¬ 
ducted  until  the  approximate  resonant  frequency  of  the 
system  was  found.  (For  selected  cases  an  exact 
determination  of  the  system  resonant  frequency  was 
subsequently  made  by  tracing,  on  an  X-Y  plotter^  the 
d-c  value  proportional  to  the  resonant  mass  accelera¬ 
tion  amplitude  as  a  function  of  the  base  excitation 
frequency.)  It  may  be  noted  here  that,  for  a  base- 
excited  single-degree  of  freedom  spring-damper-mass 
system,  resonance  occurs  at  an  angle  smaller  than  90 
degrees.  The  deviation  from  90  degrees  is  essentially 
determined  by  the  amount  of  damping  in  the  system. 

While  the  predetermined  vibrational  power  dissipation 
level  in  the  elastomer  test  sample  was  maintained  by 
adjustment  of  the  shaker  power  input  level,  the  vibra¬ 
tion  frequency  was  reduced  until  the  phase  angle  de¬ 
creased  to  a  value  between  15  and  25  degrees . 

Provided  none  of  the  acceleration  and  displace¬ 
ment  signals  showed  signs  of  abnormalities 
(distortions,  or  indications  of  nonaxial  motion  of  the 
resonant  mass),  a  data  point  was  recorded.  The  recording 
of  one  data  point  consisted  of  one  complete  printout  of 
all  amplitudes  (up  to  six  individual  signals),  associated 
phase  angles  and  test  frequency,  and  all  temperatures 
from  the  thermocouples  located  in  the  test  specimen 
holding  plates. 

Stepwise  increases  in  vibration  frequency  were  imposed 
and  between  six  and  eight  data  points  were  recorded 
until  the  phase  angle  reached  approximately  165  degrees. 

Tests,  comprising  steps  1  through  6,  were  then  repeated 
with  each  of  the  remaining  masses  in  turn,  each  mass 
giving  a  dynamic  system  with  a  different  resonant 
frequency. 

Shake  table  power  was  turned  off  and,  for  the  compression 
tests  the  preload  on  the  elastomer  sample  was  increased 
to  the  next  higher  value. 

Steps  4  through  7  were  repeated  until  data  at  all  desired 
preloads  were  obtained. 

Shaker  table  power  was  turned  off  and,  for  the  compression 
sample,  all  preload  on  the  elastomer  removed.  The  elastomer 
test  sample  was  then  given  time  (at  least  one  hour)  to 
recover  its  original,  uncompressed  height  (as  indicated 


by  displacement  probe  readings)  before  tests  at  either 
higher  ambient  temperatures  or  a  higher  power  dissipa¬ 
tion  level  were  commenced. 

The  result  of  this  test  sequence  is  a  printed  data  set  defining  the 
measured  amplitudes  for  each  vibration  sensor,  and  the  temperature 
values  at  each  frequency  point  for  a  series  of  mass  values.  Together 
with  data  specifying  the  test  conditions,  this  data  set  is  entered 
directly  onto  punched  cards  which  are  processed  by  a  computer.  All 
required  editing  and  averaging  functions  are  performed  by  the  computer 
then  stiffness  and  damping,  K  ,  are  calculated  at  each  frequency 
as  follows: 


2  2 

CO  M  [cy  •^0'  cos  cp] 

(66) 

2 

a  2a  cos  cp  +  1 

2  „ 

CO  M  sin  cp 

(67) 

2 

a  ^  2a  cos  cp  +  1 

where  a  is  the  transmissibility ,  that  is  the  ratio  of  absolute 
values  of  output  to  input  acceleration;  and 
Cp  is  the  phase  angle  between  input  and  output  acceleration 
signa Is . 

The  development  of  these  expressions  for  stiffness  and  damping  is  pre¬ 
sented  in  Appendix  B. 

Test  Parameter  Ranges 

It  is  the  purpose  of  this  Section  to  provide  an  overall  guide  to  the 
test  cases  reported  for  this  investigation  and  to  indicate  the  range 
of  parameters  covered  in  each  case.  The  method  of  presentation  chosen 
is  predominantly  that  of  data  plots  showing  calculated  elastomer  stiff 
ness  and  damping  (in  N/m)  versus  test  frequency.  These  plots  were 
originally  prepared  by  computer  line-printer,  and  offered  a  rapid, 
economical  means  to  evaluate  the  combined  effect  of  such  constant 
test  parameters  as  dissipation  power  level  in  the  specimen  and 
ambient  temperatures.  The  purpose  of  presenting  these  plots  in 
this  section  is  to  indicate  the  extent  of  parameters,  particularly 
in  the  frequency  domain.  Lines  through  individual  data  points  have, 
therefore,  not  been  extrapolated  beyond  either  the  lowest  or  the 
highest  frequency  data  point.  (The  following  section  of  the  report 
presents  the  test  data  in  expanded  format,  one  frame  for  each  test 
condition,  together  with  theoretical  correlations  and  comparisons.) 

Tests  Performed  on  Shear  Samples 

Tests  were  performed  on  two  different  shear  samples.  One  sample  con¬ 
sisted  of  eight  elastomer  strips,  each  1.27  cm  high  and  the  second 
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sample  was  made  up  of  four  strips  each  having  twice  the  height 
(2.54  cm)  of  the  eight-strip  group.  The  eight-strip  sample  was 
tested  over  a  range  of  ambient  temperatures  and  the  four-strip 
sample  was  tested  over  a  range  of  power  dissipation  levels.  Test¬ 
ing  of  the  eight-strip  sample  was  suspended  because  it  was  believed 
that  chemical  decomposition  of  the  elastomer  had  been  initiated 
when  the  ambient  temperature  had  befen  allowed  to  go  up  to  100-110  C 
during  the  93  C  tests.  Evidence  of  damage  was  inconclusive  -  a 
retest  at  32  C  (Fig,  40)  did  not  show  significantly  different  be¬ 
havior  from  original  tests,  but  subsequent  tests  at  5  C  ambient 
temperature  showed  elastomer  damping  behavior  that  appeared 
anomalous  when  compared  to  data  obtained  at  other  temperatures. 


The  range  of  test  parameters  (frequency  and  ambient  temperature) 
covered  in  tests  of  the  eight-strip  sample  are  shown  in  Figure  41. 

The  ranges  of  power  level  and  frequency  covered  in  teats  of  the 
four-strip  sample  are  shown  in  Figure  42.  A  listing  of  shear  test 
cases  is  given  in  Table  VI. 

The  power  dissipation  levels  are  nominal  values,  representing  the 
average  of  values  calculated,  after  the  fact,  from  measured  frequency, 
damping  and  amplitude  (power  =  0.5  m  Y*  K2).  The  values  for  shear 
and  compression  specimens  are  generally  higher  than  planned.  It  was 
Intended  that  each  test  would  be  conducted  at  a  constant  power  level 
was  identified,  towards  the  end  of  the  series  of  shear  speci¬ 
men  tests,  that  a  constant  power  meter  setting  actually  resulted  in 
a  variation  with  frequency:  low  frequency  dissipation  was  higher  than 
average  and  high  frequency  dissipation  was  lower  than  average.  The 
average  value  was  also  higher  than  Intended.  This  characteristic  is 
shown  in  Figure  43.  Following  a  recalibration  of  the  meter  it  was 
possible,  during  the  cartridge  tests,  to  maintain  more  consistent, 
targeted,  dissipation  levels- 

The  highest  test  frequency  is  limited  by  the  lightest  resonant  mass 
without  the  preload  piston.  This  limit  is  0.5  Kg  and  frequencies  for 
the  stiffest  shear  specimen  are  as  high  as  1600  Hz.  The  lowest  test 
frequency  is  governed  by  the  highest  mass  which  can  be  resonated  with¬ 
out  rocking  motion  and  without  applying  an  excessive  preload  (arbitrar¬ 
ily  limited  to  one  percent).  The  normal  upper  mass  limit  was  15  Kg. 

It  may  be  seen,  in  Figures  41  and  42  that,  with  these  limitations,  the 
test  frequency  range  covered  shifts  downwards  as  the  stiffness  of  the 
test  specimen  falls  -  either  with  increasing  temperature  or  with  in¬ 
creasing  dissipation. 

Figures  44,  45  and  46  present  some  interesting  comparisons  which  relate 
mainly  to  the  test  method.  In  Figure  44  test  series  4-01->06  and 
4-19 >24  are  compared.  The  conditions  for  each  set  were  Identical  but 
the  two  tests  were  separated  by  a  two  week  time  delay  in  which  the  rig 
was  disassembled  and  reassembled.  The  agreement  between  the  two  is 
good  and  indicates  that  good  repeatability  can  be  achieved  with  the 
test  method. 
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Fig.  42  Stiffness  and  Damping  for  Four-Strip  Elastomer  Specimen 
at  32  C  and  Six  Power  Dissipation  Levels 
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In  Figure  45  data  sets  4-38-M3  and  4-62^67  are  compared.  The  first 
set  was  obtained  without  the  air  cylinder  (preload  piston)  and  the 
second  set  was  obtained  with  the  air  cylinder  in  place,  and  sealed,  but 
not  pressurized.  In  other  words,  apart  from  the  presence  of  the 
piston,  the  conditions  are  identical.  The  differences  between  the 
two  sets  of  data  are  small,  but  they  are  undeniably  in  the  direction 
of  lower  stiffness  and  damping  with  the  cylinder  in  place.  This 
result  is  surprising,  since  analysis  of  the  test  rig  mechanics  in¬ 
dicates  that  the  air  spring  stiffness  and  damping  will  both  tend  to 
act  in  parallel  with  the  test  specimen.  However,  the  result  was 
fully  confirmed  by  data  sets  04-01-^-06  or  04-19^24  and  04-31->35  (not 
shown  in  graphical  presentation).  This  anamolous  trend  is  not  fully 
understood  at  this  time,  but  is  assumed  to  be  the  result  of  flexibil¬ 
ity  in  test  rig  elements  (for  example,  the  piston  or  cylinder  head). 

In  Figure  46  data  sets  4-38^43  and  4-44-4-49  for  32  C  are  compared  and 
data  sets  4-50^55  and  4-56^61  for  66  C  are  compared.  In  each  com¬ 
parison  the  preload  piston  is  present,  but  between  the  first  and  second 
set,  the  only  difference  is  that  vent  holes  in  the  cylinder  have  been 
opened.  The  effect  of  this  change  is  clearly  negligible,  leading  to 
the  conclusion  that  anomalous  behavior  of  the  air  cylinder  is  not  a 
result  of  pressure  fluctuations. 

These  anomalies  with  the  preload  piston  confirm  that  small  uncertainties 
are  introduced  in  tests  with  the  piston  and  that  the  most  reliable 
elastomer  specimen  data  will  be  obtained  without  the  preload  piston. 

Tests  Performed  on  Compression  Samples 

Under  Phase  III  all  compression  sample  tests  were  performed  on  an 
assembly  of  ten  elastomer  specimens  12.7  mm  in  diameter  and  6.35  mm  in 
height  (see  Fig.  35) .  Tests  were  generally  conducted  at  2-1/2  percent 
preload  and  then  repeated  (for  each  series  of  masses)  at  five  percent 
preload.  At  a  nominal  power  dissipation  level  of  0.062  watt/cm^,  the 
compression  sample  was  tested  at  three  ambient  temperatures;  2,  66  and 
93  C.  At  32  C  ambient  temperature  the  sample  was  tested  at  increased 
power  dissipation  levels  of  approximately  0.50,  1.1  and  2.2  watt/cm^. 
(Results  from  tests  at  32  C  and  low  nominal  power  dissipation  level 
(0.062  watt/cm  )  were  available  from  a  previously  conducted  test  series 
in  Phase  II.) 


The  frequency  range  covered  in  compression  sample  tests  extended  from 
approximately  150  Hz  for  high  temperature  tests  to  1000  Hz  for  low 
temperature  tests.  Coverage  at  the  high  frequency  end  was  generally 
not  as  high  as  for  the  shear  sample  tests,  due  to  the  preload  piston 
which  Increased  minimum  mass, 

lest  range  parameters  for  all  compression  sample  teats  are  shown  In 
Table  VII. 


Graphical  representations  of  the  test  ranges  covered  for  the  compression 


TABLE  VII.  LISTING  OF  COMPRESSION  TEST  SAMPLE  CASES 


Test  Case 
Number 

Mass  On  Elastomer 

(Kg) 

Ambient 

Temper¬ 

ature 

(C) 

Preload 
On  Sample 

a) 

Power 

Dissipation 

Level 

3 

(watt /cm  ) 

10-01-05 

0.91,  1.81,  2.78,  3.73,  5.35 

66 

2.5 

0.062 

10-06-09 

0.91,  1.81,  2.78,  3.73,  5.35 

60 

5 

G.062 

10-10-13 

0.91,  1.81,  2.78,  3.73,  5.35 

93 

2.5 

O.OO,' 

10-14-17 

0.91,  1.81,  2.78,  3.73,  5.35 

93 

5 

0.062 

10-18-21 

5.35,  3.72,  1,81,  0.91 

32 

2.5 

1 

0.44 

10-32-36 

0.41,  1.81,  3.73,  5.35 

32 

5 

0,55 

10-25-28 

0.41,  1.81,  3.73,  5.35 

32 

2.5 

l.I 

10-29-32 

0.41,  1.81,  3.73,  5.35 

32 

5 

1.1 

10-37-39 

0.41,  1.81,  3.73,  5.35 

32 

2.5 

2.2 

10-41-44 

0.41,  1.81,  3.73,  5.35 

2 

2.5 

0.062 

10-45-48 

0.41,  1.81,  3.73,  5.35 

2 

5 

0.062 

sample  are  shown  in  Figures  47  and  48.  Both  figures  are  for  the  2-1/2 
percent  preload  test  case  only;  very  similar  results  were  obtained 
with  five  percent  preload  on  the  test  sample.  The  consistency  of  the 
stiffness  and  damping  values  obtained  for  the  compression  test  sample 
appears  to  be  good,  following  very  much  the  trend  lines  already  observed 
and  discussed  for  the  shear  test  sample. 

Test  Performed  on  Cartridge  Samples 

Two  cartridge  type  samples  were  tested,  one  with  elastomer  elements 
9.55  mm  wide  and  a  second  with  elements  half  as  wide  (4.75  mm).  Both 
samples  were  first  tested  at  four  ambient  temperatures:  32,  43,  66  and 
80  C,  and  then  at  several  power  dissipation  levels  with  ambient  tempera¬ 
tures  controlled  at  32  C.  The  frequency  range  covered  in  these  tests 
is  indicated,  for  each  ambient  temperature  and  power  dissipation  level 
test  case,  in  Figures  49  and  50  for  the  wide  cartridge;  and  in  Figures 
52a  &  52b  for  the  narrow  cartridge.  Due  to  the  lower  overall  stiffness 
of  the  narrower  element  sample  the  test  frequency  range  did  shift 
somewhat  towards  lower  frequencies. 

The  effect  of  static  preload  produced  by  the  air  cylinder  (2-1/2  and  5 
percent  deflection  of  the  height  of  the  elastomer  material  annulus)  was 
investigated  for  the  wide  cartridge  specimen  only.  In  Figure  51  results 
without  the  preload  piston  attached  are  compared  with  results  obtained 
with  the  piston  attached  and  pressurized  to  both  2-1/2  and  5  percent 
preload.  The  effect  of  introducing  the  piston  appears  to  be  much  stronge 
than  the  effect  of  doubling  the  preload.  While  no  direct  comparison, 
showing  the  effect  of  the  piston  alone,  is  available  for  the  cartridge 
specimen,  the  indications  are,  just  as  for  the  shear  specimens,  that 
the  introduction  of  the  piston  reduces  the  observed  stiffness  and  damping 

A  listing  of  all  cartridge  sample  test  cases  is  given  in  Table  VIII. 

Static  Tests 

Static  stiffness  tests  were  conducted  with  the  elastomer  samples 
mounted  in  the  test  rig  and  the  displacement  probe  output  used  for 
^^bomatic  recording  of  elastomer  specimen  deflection.  Using  the  same 
environmental  control  as  for  dynamic  tests,  static  deflection  data 
was  obtained  at  temperatures  between  2  and  93  C.  Application  of  load 
was  controlled  manually  by  pressurization  of  the  preload  cylinder. 

Static  deflection  tests  of  elastomer  materials  differ  drastically 
from  those  associated  with  most  engineering  materials.  An  elastomer 
under  load  will  deflect  rapidly  to  a  certain  extent  as  the  load  is 
applied  and  then  continue  deflecting  at  a  much  lower  rate  for  an 
appreciable  amount  of  time.  The  rate  of  deflection  decreases  with 
time  after  each  step  in  load  application,  but  for  the  presently  con¬ 
sidered  material,  increases  with  Increased  total  load  on  the  sample. 

One  example  of  a  typical  load  deflection  relationship,  as  obtained  for 
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Fig.  49  Stiffness  and  Danqplng  for  Wide  (9.55  nm)  Cartridge  lement 
Specimen  at  Ambient  Temperatures  of  32,  43,  66  and  30  C 
and  0.044  watt/cm^  Power  Dissipation 
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Fig.  52b  Stiffness  and  Damping  for  Narrow  (4.75  mm)  Eiemeat 

Cartridge  at  32  C  Ambient  Temperature  and  5  _  e'^\:ent  Power 
Dissipation  Levels 


TABLE  VIII.  LISTING  OF  CARTRIDGE  TEST  SAMPLE  CASES 


Test  Case 
Number 

Ambient 

Temper¬ 

ature 

(C) 

4C- 01-06 

32 

4C-07-13 

43 

4C- 14-20 

66 

4C-21-26 

32 

4C-27-32 

32 

4C-33-38 

32 

4C-39-44 

32 

4C- 45-50 

80 

*4C-51-54 

32 

**4C-55-58 

32 

Mass  On  Elastomer 
(Kg) 

Sample  Element  Width  9.55  mm 
I  2.06,  2.98,  3.96,  6.51,  11.27,  15.65 


4C2 -01-04 

32 

4C2- 05-08 

43 

4C2-09-12 

66 

4C2-13-16 

80 

4C2-17-20 

32 

4C2-21-24 

32 

4C2-25-28 

32 

4C2-29-32 

32 

4C2-33-36 

32 

2.26,  3.46,  4.46,  11.77 
2.56,  4.46,  11.77,  16.14 

Sample  Element  Width  4.75  mm 

2.06,  3.96,  6.51,  15.67 


*With  Piston  (2-1/2%  Preload). 
**With  Piston  (5%  Preload). 


Power 

Dissipation 

Level 

(watt/cm^) 


0.044 

0.044 

0.044 

0.16 

0.5 

1.0 

1.8 

0.044 

0.044 

0.044 


0.022 

0.022 

0.022 

0.022 

0.11 

0.22 

0.44 

0.88 

1.32 


the  elastomer  cartridge  sample  at  66  C,  is  shown  in  Figure  53.  Beginning 
at  time  zero  with  no  load  on  the  sample  (upper  left-hand  corner  in  Figure 
53)  the  output  of  the  displacement  sensor  (digital  voltmeter  reading)  is 
high,  indicating  a  large  gap  between  sensor  and  target.  As  the  air  pres¬ 
sure  in  the  loading  cylinder  on  the  test  rig  is  increased,  first  to  4 
psi  and  then  successively  in  4  psi  steps  to  16  psi,  the  voltmeter  read- 
ing  decreases  sharply  at  first  and  then  much  more  slowly.  Approximately 
twelve  minutes  after  application  of  each  load  step  the  plotter  pen  was 
brought  back  to  the  left  side  of  the  paper  and  the  pressure  in  the  load¬ 
ing  cylinder  was  increased. 

From  Figure  53  and  similar  data  curves  the  plots  for  the  shear  ar<d  cartridg 
samples  -  Figures  54  and  55  -  have  been  assembled.  The  load  was  pi.o,■^ed 
against  the  last  deflection  value  for  each  load.  Similar  loading  was  cp- 
plied  to  the  compression  specimen  but,  in  this  case,  a  mechanical  diai. 
indicator  was  used  for  deflection  measurement.  Figure  56  shows  compression 
specimen  deflection  data. 

Examination  of  the  load  versus  deflection  curves  for  the  shear  sample 
(Fig.  54)  shows  several  trends  which  may  also  be  observed  for  the  com¬ 
pression  sample  tests  (Fig.  56),  but  are  less  pronounced  for  the  car¬ 
tridge  sample  (Fig.  55). 

1.  The  load  versus  delfection  curve  is  nonlinear  at  small 
deflections  only. 

2.  The  degree  of  nonlinearity  decreases  with  decreasing 
elastomer  stiffness. 

3.  Elastomer  stiffness  decreases  with  increasing  ambient 
temperature . 

Details  of  the  nonlinear  portion  of  the  stiffness  curve  for  the  shear 
sample  (Pig.  54)  were  confirmed,  with  greater  precision,  through  sub¬ 
stitution  of  deadweight  loading  for  air  cylinder  pressurization. 

For^the  cartridge  sample  (Fig.  55)  the  deflection  curves  do  not  seem 
to  straighten  out'^  within  the  test  range  as  they  do  for  the  other  sanqjles. 
The  limited  range  of  cartridge  sample  deflections  (about  0.33  mm),  may 
be  directly  responsible,  which  was  the  result  of  displacement  sensor 
limitations. 
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^  FIr.  56  Static  Deflection  vs  Force  for  Compression  Element  Specimen 

at  Various  Temperatures 


TEST  RESULTS 


The  teste  described  in  the  previous  sections  generated  discrete  values 
for  stiffness  and  damping  of  shear,  cylindrical  and  ring  cartridge 
elastomer  elements  at  different  values  of  frequency,  temperature  and 
dissipation  level.  The  purpose  of  this  section  is  to  present  this 
data  in  its  entirety,  as  well  as  to  provide  an  appropriate  summary  of 
it.  The  data  presented  forms  the  basis  for  further  interpretation. 

The  basic  medium  of  presentation  is  the  plot  of  stiffness  and  damping 
versus  frequency.  One  frame  is  presented  for  each  set  of  conditions 
of: 

•  Geometry 

•  Temperature 

•  Dissipation  Level 

•  Preload  (where  appropriate) 

On  each  frame  discrete  stiffness  and  damping  values  are  plotted,  as 
well  as  a  line  representing  the  least  squares  best  fit  power  law  re¬ 
lationship  between  stiffness,  damping  and  frequency,  in  the  form: 

CD  (68) 

®2 

CD  ^  (69) 

The  shear  specimen  power  law  expressions  for  different  temperatures 
are  used  to  develop  a  catalog  of  material  properties  for  this  material. 
The  material  properties  catalog  is  subsequently  used  as  part  of  the 
basis  for  predicting  stiffness  and  damping  of  the  cylindrical  compres¬ 
sion  and  cartridge  specimens. 

A  brief  description  of  the  various  sets  of  test  results  follows. 

Shear  Specimen  Tests  at  Different  Temperatures 

These  tests  were  performed  with  the  eight-element  shear  specimens, 
each  element  having  a  dimension  of  1.27  cm  in  the  direction  of  loading. 
The  tests  were  performed  with  the  simplest  configuration  possible  -  a 
mass  symmetrically  supported  by  the  shear  elements  which  are  mounted 
on  the  table.  The  preload  piston  which  has  been  used  in  previous 
tests  was  not  used,  eliminating  any  effects  of  tare  damping. 

The  data  is  presented  in  Figures  57  through  60  for  temperatures  of  32, 
50,  66  and  93  C.  The  consistent  effect  of  temperature  is  to  reduce 
stiffness  and  damping.  The  data  shows  reasonable  consistency  for  the 
three  higher  temperatures,  but  at  32  C  some  scatter  is  encountered  at 
the  highest  frequency.  The  power  law  relationship  provides,  in  each 
case,  a  reasonable  fit  to  the  data  and,  certainly,  no  more  sophisticated 
representation  would  be  justified. 


106 


107 


ig.  58  Shear  Specimen  at  Zero  Preload,  50  C  (Eight^ElCTients  as 
Shown)*  Nominal  Dissipation  0,062  watts/cm^ 
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Material  Properties  Catalog 

Applying  to  the  shear  specimen  test  data  the  interpretive  formula: 


h 

nA 


(70) 


nA 


(71) 


where  h  is  the  strained  dimension  (0.3175  cm), 
nA  is  the  total  stressed  area  (49.55  cm^). 

A  power  law  relationship  between  storage  and  loss  moduli  (G%  G and 
frequency  has  been  obtained,  as  shown  graphically  in  Figures  61  ano  •<> 
The  coefficient  values  are  shown  in  Table  IX. 


TABLE  IX.  MATERIAL  PROPERTIES  CATALOG  FOR 
POLYBUTADIENE,  NICHOLS  NEX  156G 


Temperature 

Storage  Modulus 

Loss  Modulus 

32  C 

43  C 

66  C 

93  C 

G'  =  3.686  X  10^ 

G'  =  1.904  X  10^  cu°-2627 

G'  =  1.097  X  10^ 

G'  =  1.197  X  10^ 

G"  =  8.333  X  10^ 

G"  =  3.971  X  10^  ^-0.06372 

G"  =  2.423  X  10^  ^“0.03011 

r"  ^  nc/  +0.2246 

G  =  1 .064  X  10  CD 

Together  with  the  geometry  factors  developed  in  preceding  sections, 
this  material  properties  data  has  been  applied  in  the  prediction  of 
stiffness  and  damping  of  the  cylindrical  and  cartridge  specimens.  For 
simplicity  of  ajjplication,  it  is  convenient  to  have  available  the 
values  of  G  ,  G  at  100  and  1000  Hz,  and  these  values,  based  on  Table 
IX  are  presented  in  Table  X. 


TABLE  X.  VALUES  OF  STORAGE  AND  LOSS  MODULI  FOR 
POLYBUTADIENE,  NICHOLS  NEX  156G  AT 
100  AND  1000  HERTZ 


Temperature 

G'  100  Hz 
N/m2 

G'  1000  Hz 
N/m^ 

g"  100  Hz 
N/m^ 

G"  1000  Hz 
N/m2 

32  C 

43  C 

66  C 

93  C 

1.369  X  10^ 
1.035 

0.7755 

0.5803  vl' 

2.189  X  10^ 
1.894 

1.5606 

1.0201  i 

3.659  X  10^ 
2.633 

1.996 

1.2297  vl/ 

! 

2.728  X  1 

2.274 
1.862 
1.2949  ^ 
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Fig,  62  Influence  of  Temperature  on  Loss  Modulus  (Shear  Spec 
Data»  Eight  Elements) 


Compression  Specimen  Tests  at  Different  Temperatures 

This  series  of  tests  was  performed  with  the  ten-element  compression 
specimens,  each  element  1.27  cm  in  diameter  and  0.635  cm  high.  In 
this  case,  the  preload  piston  was  a  necessary  feature  of  the  tests  to 
satisfy  the  requirement  that  the  compression  elements  not  bs  allowed 
to  experience  tension.  Geometrical  preload  values  of  2.5  and  5  per¬ 
cent  were  applied.  Temperature  values  tested  were  66,  93  and  5  C. 

These  test  results  are  presented  in  Figures  63  through  68.  Again, 
reasonable  consistency  between  the  test  data  and  the  fitted  line  was 
obtained.  The  consistent  trend,  as  for  the  shear  specimens,  was  that 
stiffness  and  damping  fall  with  increasing  temperature. 


Included  in  the  figures  for  66  and  93  C  are  predictions  of  stiffness 
as  a  function  of  frequency  using  the  empirical  geometry  coefficients 
developed  on  the  basis  of  previous  (Phase  IT)  test  data  and  using 
the  material  property  values  obtained  from  the  present  shear  specimen 
tests,  as  presented  in  Table  IX.  In  ajidition,  predictions  using  the 
static  shape  coefficient  value  (3  =3  =2)  are  presented  on  this 

figures . 

Very  similar  predictions  are  obtained  with  either  approach.  The  pre¬ 
dictions  of  stiffness  are  within  ten  percent  of  measurement  at  66  C  and 
within  twenty  percent  at  93  C.  Predictions  of  damping  are  within  25  per¬ 
cent  at  66  C  and  within  twenty  percent  at  93  C.  No  predictions  are  made 
at  5  C  since  material  property  data  was  not  generated  at  this  temperature. 

While  the  stiffness  predictions  show  similar  trends  to  the  measurements, 
in  terms  of  rate  of  increase  with  frequency,  the  damping  predictions  in¬ 
dicate  a  slight  decrease  with  frequency  and  the  measurements  reveal  a 
more  pronounced  Increase  with  frequency. 

Compression  Specimen  Tests  at  Different  Dissipation  Rates 

These  tests  were  performed  with  the  ten-element  compression  specimens 
at  32  C.  Preloads  of  both  2,5  and  5  percent  were  investigated  by  using 
the  preload  piston. 

The  variation  of  stiffness  and  damping  with  frequency  at  different 
nominal  dissipation  levels  is  presented  in  Figures  69  through  73.  The 
corresponding  variation  of  amplitude  is  shown  in  Figures  74  and  75  and, 
using  these  two  sets  of  data,  the  details  of  each  test  point  in  terms 
of  amplitude,  frequency  and  temperature  may  be  recreated. 

The  general  effect  of  increasing  dissipation  is  to  decrease  both  stiff¬ 
ness  and  damping,  at  all  frequencies,  for  these  compression  specimens. 

As  previously  discussed,  it  is  hypothesized  that  a  major  reason  for 
this  degradation  is  the  increase  in  the  material  temperature.  The 
comparison  between  this  measured  data  and  predictions  of  the  simple 
thermoviscoelastic  model  will  be  presented  in  the  next  Section  of  the 
report. 
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Fig,  63  Stiffness  and  Damping  vs  Frequency,  Compression  Specimen 
(Ten  Elements  with  Dimensions  as  Shown).  2-1/2  Percent 
Preload,  66  C.  Nominal  Dissipation  0.062  watts/cm^. 
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Stiffness  and  Damping  vs  Frequency.  Compression  Specimen. 
Five  Percent  Preload,  66  C.  (Ten  Elements  with  Dimensions 
as  Shown).  Nominal  Dissipation  0.062  watts/cm^.  2-1/2 
Percent  Preload. 
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Hg.  65  Stiffness  and  Damping  ■vs  Frequency#  Compression  Specimen# 
2-1/2  Percent  Preload,  93  C.  Nominal  Dissipation  0.062 
watts/cm'* 


STIFFNESS  AND  DAMPING,  N/m 


FREQUENCY,  H* 


Fig,  67  Stiffness  and  Damping  vs  Frequency.  Compression  Specimen. 
5  C.  2^1/2  Percent  Preload;  Nominal  Dissipation  0.062 
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Fig,  68  Stiffness  and  Damping  vs .Frequency.  Compression  Specimen. 

5  C,  5  Percent  Preload;  Nominal  Dissipation  0.062  watt/cm^. 
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Fig.  69  Stiffness  snd  Damping  vs  Frequency*  Compression  Specimen* 
32  C.  Nominal  Dissipation  0.4A  watt/cm^ 
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Fig.  70 


Damping  vs  Frequency.  Compression  Specimen. 
2-1/2  Percent  Preload,  32  C.  Nominal  Dissipation  1.1  watts/ci. 
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Fig.  71  High  Dissipation  Stiffness  and  Damping  vs  Frequency* 
Compression  Specimen.  2-1/2  Percent  Preload,  32  C. 
Nominal  DJeslpatlon  2.2  watts/cm^ 
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Fig.  73  Stiffness  and  Damping  vs  Frequency.  Compression  Specimen. 

5  Percent  Preload,  32  C.  Nominal  Dissipation  1.1  watts/cm^ 
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Amplitude  as  a  Function  of  Frequency.  Compression  Specimen. 

5  Percent  Preload,  32  C  (Ten  Elements  with  Dimensions  as  Shown) 


Shear  Specimen  Tests  For  Different  Dissipation  Rates 


These  tests  were  performed  with  four-element  shear  specimens.  Each 
element  is  2.54  cm  long  in  the  direction  of  loading,  the  strained  di¬ 
mension  is  0.3175  cm  thick,  and  the  width  of  the  elements  is  4.877  cm. 
As  with  the  eight-element  shear  specimen,  it  is  possible  to  perform 
these  tests  without  the  uncertainties  introduced  by  the  presence  of  the 
preload  piston.  All  the  high  dissipation  tests  were  performed  at  a 
temperature  of  32  C. 

As  a  reference,  a  low  dissipation  level  of  0.035  watt/cm^  has  been 
tested,  followed  by  nominal  dissipation  values  5,  10,  20,  50  and  125 
times  the  reference.  The  data,  in  terms  of  stiffness  and  damping,  is 
presented  in  Figures  76  through  81  and  the  amplitudes  as  a  function  of 
frequency  in  Figure  82. 


Reasonable  consistency  between  the  test  data  and  the  fitted  line  is  ob¬ 
tained,  with  scatter  of  the  order  of  ten  percent  of  the  stiffness  and 
damping  level.  Scatter  increases  at  the  highest  dissipation  levels. 

It  is  generally  a  trend  that  both  stiffness  and  damping  fall  with  in¬ 
creasing  dissipation  level. 

The  slope  of  the  lines  also  tend  to  increase  with  increasing  dissipation 
levels  -  indicating  a  more  pronounced  effect  of  high  dissipation  at  low 
requencies  .  This  could  result  from  the  power  meter -characteristic  dis- 
cussed  in  the  previous  Section. 

Cartridge Tests  at  Different  Temperatures 

These  tests  were  preformed  on  two  different  cartridge  test  specimens 
^ch  consisted  of  four  elements  mounted  in  the  same  fixtures  and  holders. 

e  inner  and  outer  diameters  of  the  elements  were  the  same  for  the  two 
specimens  (1.90  cm  and  2.86  cm);  the  difference  was  the  length  -  0  47  cm 
in  one  case  and  0.95  cm  in  the  other.  Data  for  the  shorter  elements  are 
presente  in  Figures  83  through  86  and  for  the  larger  elements  in  Figures 
»/  through  90.  The  temperature  range  covered  is  32  to  80  C.  The  enerev 
dissipation  rate  is  0.022  watts/cm^. 

generally  consistent  with  the  fitted  power  law  line,  except 
that  at  low  frequencies,  for  the  shorter  length  elements,  the  damping 
values  fall  with  frequency  before  starting  to  increase  with  frequency 
above  150  to  200  Hz.  The  influence  of  increasing  temperature  is  to  de- 
crease  the  stiffness  and  damping.  In  general,  the  influence  of  frequency 
is  to  increase  both  stiffness  and  damping  at  a  similar  rate. 

mefboH^nr'^  on  these  test  data  plots  are  predictions  using  the  various 
methods  of  prediction  for  the  cartridge  geometry  previously  discussed. 

In  all  cases,  the  Gobel  predictions  and  Beam-Column  predictions  are  in¬ 
cluded.  In  most  cases  the  plane  stress  line  is  also  included. 


6x10^ 


Flp,,  76  Stiffness  and  Damping  vs  Frequency.  Shear  ioecimens 

Zero  Preload,  32  C.  Nominal  Dissipation  0.035  watts/cm^ 
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Fig.  78  Stiffness  and  Damping  vs  Frequency.  Shei.  Sfecimens 
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;  vs  i<requency.  Shea-  Spicimens. 
Nominal  Dissipatica  1.76  watts/cm 


Amplitude  as  a  Function  of  Frequency.  Shear 
Zero  Preload,  32  C 
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Fig.  84  Stiffness  and  Damping  for  Cartridge  Specimen. 
Nominal  Dissipation  0.022  watts/cm^ 
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Fig.  86  Stiffness  and  Damping  for  Cartridge  Specimen 
Nominal  Dissipation  0.022  watts/cm^ 


KIr.  87  Stiffness  and  Damping  for  Cartridge  Specimen 
Nominal  Dissipation  0.022  watta/cm^ 
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Dynamic  Stiffness  and  Damping  for  Cartridge  Sp;  Imen.  44  c. 
Nominal  Dissipation  0.022  watts/cm  (Four  FT  aments  as  Shown) 
(D^ensions  as  Shown)  . 


ELASTOMER 


Damping  for  Cartridge  Specimen.  80 
022  watts/cm^  (Four  Elements  as  Sho' 


With  some  exceptions,  for  damping  of  the  shorter  specimen,  the  test 
data  lies  within  the  band  defined  by  the  highest  and  lowest  predic¬ 
tion  lines.  The  shorter  specimen  data  generally  approaches  the 
Beam-Column  predictions  most  closely,  while  the  longer  specimen  data 
more  generally  lies  midway  between  the  Beam-Column  and  Gobel  predic- 
tion  lines. 


The  stiffness  data  shows  a  similar  slope  to  the  predictions,  whereas 
the  damping  test  data  shows  a  different  trend  from  the  predictions. 

Ihe  material  properties  for  the  elastomers  obtained  from  shear  speci¬ 
men  tests  indicate  a  slight  decrease  in  loss  modulus  with  frequency. 
Since  the  cartridge  damping  predictions  are  related  to  the  loss  modu¬ 
lus  via  a  frequency  independent  geometry  factor  they  also  reflect  this 
decrease  with  frequency.  The  cartridge  test  data,  however,  shows  a 
consistent  incresse  with  frecjuency. 


Cartridge  Tests  at  Different  Dissipation  Levels 

These  tests  are  performed  with  the  same  two  cartridge  test  specimens 
as  for  the  tests  at  different  temperatures.  The  temperature  for  these 
tests  was  32  C.  Data  for  the  shorter  specimen  are  presented  in  Figures 
91  through  96  and  for  the  larger  specimen  in  Figures  97  through  101. 

In  Figures  96  and  101,  corresponding  to  each  plot  of  stiffness  and  damp¬ 
ing  versus  frequency,  is  a  line  showing  amplitude  as  a  function  of 
frequency.^  With  one  exception,  there  is  good  distinction  between  the 
lines  on  Figures  96  and  101.  However,  in  Figure  96,  above  210  Hz 
the  points  identified  as  1.32  watts/cm^  nominal  dissipation  lie  very 
close  to  the  line  for  0.88  watts/cm^.  Apparently  the  required  test 
conditions  for  1.32  watts/cm  were  not  maintained  above  210  Hz.  Com¬ 
parison  ol  Figures  94  and  95  confirm  that,  below  210  Hz,  there  is  sig- 
niticant  difference  in  the  stiffness  and  damping  data,  whereas,  above 
lix,  the  two  figures  present  very  similar  results.  Figure  95  is 
intended,  for  completeness  since,  toghether  with  the  data  of  Figure 
96  It  fully  defines  the  test  conditions  for  individual  points. 

The  effects  of  dissipation  level  are,  generally,  to  decrease  both  stiff¬ 
ness  and  damping  although  this  trend  is  less  pronounced  at  high  fre- 
quency  than  at  low  frequency.  In  other  words,  the  lines  of  stiffness 
and  damping  versus  frequency  becomes  steeper  at  high  dissipation  levels. 

Specimen  Catalog; 

The  power  law  curve  fit  coefficients  for  all  the  specimens  whose  stiff¬ 
ness  and  damping  are  plotted  in  Figures  57  through  101  are  summarized 

dissipation  and  geometry  in  Tables  XI, 

XII  and  XIII.  These  tables  form  the  Specimen  Catalog  for  shear  com¬ 
pression  and  cartridge  geometries.  ’ 


Transient  Behavior 


In  high  dissipation  tests  with  the 
dissipation  was  readied  at  which  it 
a  subre.sonant  condition.  Wlien  the 


compression  specimens,  a  level  of 
was  no  longer  possible  to  maintain 
input  acceleration  to  the  shake  table 
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Fig.  91  Stiffness  and  Damping  for  Cartridge  Specimen 
Nominal  Dissipation  0.11  watts/cm^ 


al  Dissipation  0.88  watts/ 


Amp  i  i  1  ikU'  vs  Fr ('qinnicy  .  Narrow  Cartridge  Specimen,  High 
Dissipation  Tests.  32  C  (Four  Elements  with  Detail  as  Shown) 


200  300  400  600  8001000  2000  <?'XX)  6000  10,000 

FREQUENCY,  HZ 
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Fig.  98  Dynamic  Stiffness  and  Damping  for  Cartridge  Specimen.  32  C 
Nominal  Dissipation  0.50  watts/cm^  (Four  Elements  as  Shown) 


Fig.  100  Dynamic  Stiffness  and  Damping  for  Cartridge  Specimen.  3 

Nominal  Dissipation  1.8  watts/cm^  (Four  Elements  as  Shown) 


TABLE  XI  .  SPECIMEN  CATALOG 
(SHEAR  SPECIMEN) 
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REPRODUCIBILITY  OP  THE 
ORIGINAL  PAGE  IS  POOR 


^Dimensions  are  for  thickness  (h)  ,  wicith  (b) ,  and  length  (x)  . 


TABLE  XII.  SPECIMEN  CATALOG 

(COMPRESSION  ELEMENTS) 


L 

'im 

im 


m 


Im 


157 


^Dimensions  are  diameter  (d),  and  height  (h). 


TABLE  XIII,  SPECIMEN  CATALOG 

(CARTRIDGE  ELEMENTS) 


Dimensions  are  inner  diameter  (D.),  outer  diameter  (D  ),  and  width  (i,) 


reached  a  critical  level,  a  slow  but  consistent  change  in  system 
conditions  would  set  in.  With  no  further  adjustments  to  frequency 
of  input  acceleration,  the  output  acceleration  would  steadily  increase 
with  time,  reach  a  peak,  then  reduce  and  settle  at  a  steady  value. 
During  this  change  in  acceleration  the  phase  angle  would  increase 
from  as  low  as  30^  and  stabilize  at  a  value  above  90°.  This  change 
in  phase  angle  was  a  clear  indication  that  the  operating  condition  of 
the  system  had  changed  from  subresonant  to  superresonant.  This  tran¬ 
sient  phenomenon  was  also  accompanied  by  a  significant  increase  by  as 
much  as  10  C  in  the  temperature  measured  near  the  specimen. 

With  a  higher  input  acceleration,  the  same  phenonmenon  would  ^.^cur 
but  its  time  scale  would  be  reduced.  When  just  above  the  critical 
acceleration  level,  the  peak  amplitude  would  occur  after  several 
utes.  However,  levels  of  acceleration  could  be  reached  where  the  peak 
amplitude  would  occur  within  15  seconds. 

The  phenomenon  was  not  observed  in  the  shear  specimen  tests  or  in  the 
cartridge  specimen  tests. 

Some  typical  traces  of  output  acceleration  as  a  function  of  time  are 
shown  in  Figure  102.  The  frequency  in  all  cases  is  433  Hz  and  the 
supported  mass  is  0.908  kg.  The  lowest  line  on  the  figure  illustrates 
the  result  of  an  input  acceleration  below  the  critical  level  at  which 
the  trnasient  phenomenon  sets  in.  In  the  center  line  on  this  figure 
the  input  acceleration  is  just  above  the  critical  level.  Over  several 
minutes  the  output  acceleration  increases  at  a  slow  steady  rate  and 
then  increases  sharply  until  a  peak  is  reached,  followed  by  a  gradual 
drop  to  an  asymptotic  value  slightly  lower  than  the  peak  value.  In 
the  upper  line  of  this  figure,  the  rate  of  increase  in  output  acceler¬ 
ation  is  high  from  outset.  The  peak  is  reached  within  15  seconds  but 
the  gradual  drop  in  amplitude  to  an  aymptotic  value  still  takes  over 
a  minute. 

Phase  angles  at  the  beginning  and  end  of  each  trace  are  shown  as  cir¬ 
cled  values.  The  increase  in  phase  angle  from  below  to  above  90®  for 
the  two  lines  with  a  peak  is  apparent. 

Tlie  fact  that  this  is  a  transient  phenomenon,  the  fact  that  a  chnage 
in  temperature  is  involved,  and  the  fact  that  the  time  scale  of  the 
phenomenon  is  controlled  by  the  input  acceleration  and  hence  by  the 
dissipation  rate,  strongly  suggests  that  the  cuase  is  a  time-varying 
change  in  the  thermal  state  of  the  elastomer  induced  by  internal 
heating . 


ACCELERATION 


DISCUSSION 


High  Dissipation  Results 

As  shown  by  the  test  data  of  the  preceding  section,  the  influence  of 
a  high  dissipation  rate  is  generally  to  reduce  the  effective  stiffness 
and  damping  of  an  elastomer  specimen.  An  obvious  hypothesis  is  that 
internal  heating  reduces  local  material  properties,  and,  on  the  basis 
of  this  hypothesis,  the  thermoviscoelastic  model  described  earlier  was 
developed.  The  hypothesis  can  now  be  tested  by  comparing  the  predic¬ 
tions  of  this  model  with  the  test  data, 

A  parameter  required  by  the  model  is  thermal  conductivity,  and  no  va^^e 
of  conductivity  is  available  for  the  material  tested.  However,  since 
all  other  parameters  of  the  model  can  be  assigned  values,  once  material 
properties  and  geometry  factors  are  available,  the  approach  taken  has 
been  to  exercise  the  model  for  a  range  of  reasonable  values  of  con¬ 
ductivity,  and  to  determine  the  value  giving  the  best  agreement.  The 
initial  object  of  this  comparison  is  the  compression  test  data.  In 
Figures  103,  104  and  105  stiffness  and  damping  are  plotted  against 
amplitude  for  frequencies  of  175,  380  and  600  Hz,  Superimposed  on  the 
compression  specimen  test  data  are  predictions  of  the  model  for  conduc¬ 
tivity  values  of  20  x  10“^  and  6.7  x  10""^  watts/ (cm  C).  It  may  be 
seen  that  the  predictions  for  a  conductivity  of  6.7  x  10"^  match  the 
test  data  with  reasonable  success,  both  in  terms  of  trends  and  in 
absolute  values  of  stiffness  and  damping. 

The  value  of  6.7  x  10“^  for  conductivity  is  not  unreasonable  (rubber 
conductivity  is  1.49  x  10"^)  but  there  is  no  evidence,  apart  from  the 
empirical  comparison  that  this  value  is  correct  for  polybutadiene  nor, 
therefore,  is  this  comparison  absolute  evidence  of  the  validity  of  the 
model . 

The  observed  transient  change  in  amplitude  and  phase  angle  under  cer¬ 
tain  high  dissipation  conditions,  as  presented  in  the  Results  Section, 
is  qualitatively  predictable  by  the  one-dimensional  thermoviscoelastic 
analysis  described  earlier.  As  an  exercise  for  the  model,  an  attempt 
was  made  to  reproduce  the  conditions  of  Figure  102,  and  to  predict 
comparable  transient  results.  In  addition  to  conductivity  no  value  is 
available  for  specific  heat  for  this  material,  although  a  density  value 
was  obtained  from  the  manufacturer  (0.043  Ib/in,^).  The  conductivity 
and  specific  heat  values  used  were  6.7  x  10“^  watts/(cm  C)  and  0.0622 
nm/kg  C.  The  frequency  was  433  Hz  and  the  supported  mass  was  0.908  kg. 

Figure  106  presents  some  of  the  results  of  this  prediction  for  two 
levels  of  input  acceleration  (17.25  and  24  g),  which  are  similar  to 
those  applied  to  produce  the  transient  phenomenon  in  Fig.  102. 

A  thermal  transient  phenomenon  with  a  transition  through  resonance 
is  predicted  for  both  Input  accelerations,  but  the  resonant  peaks  are 
much  sharper  than  those  found  in  a  similar  time  scale  in  the  corre¬ 
sponding  test  conditions  shown  in  Figure  102.  The  predicted  internal 
peak  temperatures  are  surprisingly  high,  although  no  measurements  of 
internal  temperature  are  available  for  comparison. 
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Fig,  106  Fxedlctlon  of  Thermal  Transient  for '  Comparison  With 
Figure  102,  433  Hz,  0,908  Kg,  Supported  Mass,  Com¬ 
pression  Specimen.  Conductivity  6.7  x  10“^  watts/cm  C 


The  inability  to  completely  recreate  the  transient  conditions  is  not 
surprising,  since  very  close  simulation  of  the  mechanical  properties 
of  the  elastomer  specimen  and  its  dependence  on  temperature  and  dis¬ 
sipation  would  be  required.  However,  the  ability  qualitatively  to 
predict  the  phenomena  on  the  basis  of  a  thermoviscoelastic  model  with 
similar  levels  of  input  acceleration  is  evidence  of  the  importance  of 
thermal  effects  in  determining  the  high  dissipation  behavior  of  the 
compression  specimens , 

A  further  check  on  the  model  was  sought  by  comparison  of  its  predic¬ 
tions  with  test  data  for  the  shear  specimen.  The  same  value  of  con¬ 
ductivity  most  successful  for  the  compression  specimens  has  been 
applied  (K  =  6.7  x  10"^  watts/cm  C) ,  Figures  107  through  110  show 
typical  results  which  exhibit  significant  differences,  both  qualitative 
and  quantitative, between  tests  and  predictions.  The  predictions  give 
much  higher  stiffness  than  the  test  data.  It  turns  out  that  the  area 
for  heat  flow  is  so  large  and  the  length  of  the  heat  flow  path  so 
small  in  the  shear  specimens,  that  temperature  levels  are  much  lower 
than  in  the  compression  specimens.  Thus,  the  predictions  based  on 
internal  heating  and  a  corresponding  reduction  of  material  property 
values  show  only  minor  deviations  from  the  isothermal  values. 

It  becomes  clear  that  some  other  mechanism  is  influential  in  reducing 
stiffness  and  damping  of  the  shear  specimens.  Reference  to  available 
literature  on  the  subject  shows  that  another  mechanism  has,  indeed, 
been  recognized  though  not,  as  yet,  analyzed  or  understood.  This 
mechanism  shows  itself  as  a  softening  of  shear  elements  at  moderate 
to  high  strains.  A  transition  region  in  amplitude  has  been  identified 
as  lying  between  strains  of  0.001  and  0.005.  Below  this  strain  level 
the  stiffness  and  damping  are  independent  of  strain.  Above  this  level 
of  strain  a  fall-off  in  stiffness  occurs.  It  is  reported  by  Warnaka 
(ref,  4)  that  consistent  behavior  is  obtained  for  different  materials. 

This  softening  effect  has  not  been  previously  reported  at  these  fre¬ 
quencies  for  compression  loading  and  there  is  apparently  some  dif¬ 
ference  in  the  behavior  of  compression  and  shear  specimens  with  respect 
to  high  dissipation.  Thermal  effects  clearly  account  for  more  of  the 
stiffness  degradation  in  compression  elements  than  in  shear  elements. 

In  an  attempt  to  clarify  the  relative  importance  of  thermal  effects 
and  the  '  strain  softening '  effect,  stiffness  and  damping  data  for  the 
compression  and  shear  specimens  has  been  plotted,  separately,  against 
strain  and  against  dissipation  rate,  as  independent  variables.  The 
results  are  shown  in  Figures  111  and  112  for  shear,  113  and  114  for 
compression.  The  compression  specimen  plots  include  data  from  the 
Phase  II  series  of  tests.  Some  interesting  observations  may  be  made: 

1.  Tlic  shear  specimen  stiffness  and  damping  data  shows 
reasonably  close  correlation  with  strain. 
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Comparison  of  Predictions  and  Measurements  of  St  lf'  less  and 
Damping  Under  High  Dissipation  Conditions.  She^r  ipecimen. 
(Four  Elements) 
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Fig.  109  Comparison  of  Predictions  and  Measurements  of  St if  ness  a: 

Damping  Under  High  Dissipation  Conditions.  She«r  ‘:^eciin’^’ 
(Four  Elements) 
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Fig.  110  Comparison  of  Predictions  and  Measurements  of  Stiffness  and 
Damping  Under  High  Dissipation  Conditions.  Shear  Specimen. 
(Four  Elements) 
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.001  .01  .i 

STRAIf^,  m/m 

Fig.  Ill  Stiffness  and  Damping  as  a  Function  of  Strain.  Shear 
Specimen.  32  C  (Four  Elements) 


DISSIPATION,  WATTS /CM^ 

Fig.  112  Stiffness  and  Damping  as  a  Function  of  Dissipation. 
Specimen.  32  C  (Four  Elements) 


STRAIN  (DOUBLE  AMPLITUDE),  m/m 


Stiffness  and  Damping  as  a  Func 
Specimen.  32  C.  (Ten  Elements) 


Damping  as  a  Function  of  Strain.  Compression 


STIFFNESS  AND  DAMPING,  N/ro 


DISSIPATION,  WATTS  /  CM^ 


Fig.  lU  Stiffness  and  Damping  aa  a  Fdnctlon  of  Dlsalpatlon. 
Compression  Specimen.  32  C  (Ten  Elements) 


2,  The  shear  specimen  data  shows  poor  correlation  with 
dissipation  rate¬ 
s'  The  compression  specimen  data  shows  reasonably  close 
correlaticr  either  with  strain  or  with  dissipation 
rates . 

4.  The  data  is  generally  at  a  high  enough  strain  rate 
that  strain  independence  is  not  apparent,  (although 
more  extensive  data  would  be  required,  at  low 
strains,  to  determine  where  strain  independence  does 
occur). 

5.  Damping  values  exhibit  a  peak,  in  the  case  of  the 
shear  and  compression  specimens,  for  strains  in  the 
region  of  one  percent. 

In  summary,  these  results  justify  the  following  statements : 

Stiffness  and  damping  of  elastomer  elements  in  the  frequency  range 
100  to  1000  Hz  are  influenced  by  two  aspects  of  vibration  amplitude: 
firstly  strain,  which  has  a  ’softening  effect’  above  a  certain  transi¬ 
tion  value;  secondly,  self-heating  due  to  thermal  dissipation,  which 
causes  reduction  in  values  of  local  material  properties.  The  rela¬ 
tive  importance  of  these  two  effects  is  governed  by  the  geometry, 
type  of  loading,  and,  probably,  by  frequency  and  material.  The  strain 
softening  effect  appears  to  dominate  for  the  shear  specimens  of  poly¬ 
butadiene. 

Further  informal  ion  will  be  required  before  the  relative  contribu¬ 
tions  of  thermal  and  strain  softening  can  be  established  for  the 
compression  specimens.  The  conductivity  must  be  separately  deter¬ 
mined  before  the  one-dimensional  thermoviscoelastic  model  can  be 
tested  further.  The  low  strain  region  must  be  more  extensively 
explored  to  verify  the  existence  of  a  strain  ’transition’  region,  to 
determine  what  conditions  are  required  to  establish  strain  indepen¬ 
dence,  and  to  establish  the  influence  of  geometry  on  the  strain 
dependence.  Temperature  measurements  within  the  material  under  high 
strains  would  provide  additional  Insight. 

The  present  indications  are  that  high  amplitude  stiffness  and  damping 
follow  closely  a  unique  relationship  with  strain.  If  this  observa¬ 
tion  can  be  supported  with  further  test  data  it  offers  an  interesting 
and  simple  approach  to  accounting  for  high  amplitude  effects  in 
design . 

Damping  Performance 

The  effect  of  strain  on  damping  performance  is  further  illustrated  in 
Figure  115.  Here  the  loss  coefficient,  T),  has  been  plotted  as  a  func¬ 
tion  of  strain.  The  loss  coefficient  is  a  good  measure  of  the 
elastomer’s  ability  to  damp  vibrations  of  a  mechanical  system  and, 


STRAIN  (BASED  ON  DOUBLE  AMPLITUDE  ) 

115  Loss  Coefficient  as  a  Function  of  Strain. 


for  a  one-degree-of- freedom  system  such  as  is  used  for  the  Base  Ex¬ 
citation  Resonant  Mass  tests,  the  damping  ratio  (£  =  C/C^)  is  equal 
to  one-half  of  the  loss  coefficient.  Consider  Figure  115:  for 
strains  below  0.02  there  is  a  universal  increase  in  loss  coefficient 
for  all  specimens;  at  a  strain  of  0.02  all  specimens  show  a  loss  co¬ 
efficient  in  the  range  0.21  to  0.26  (10.5  to  13  percent  of  critical 
damping).  Above  a  strain  of  0.02,  however,  there  is  some  divergence 
in  loss  coefficient  behavior  between  the  different  specimens.  The 
shear  specimen  reveals  a  steady  increase  in  loss  coefficient  to  a 
value  of  0.3  at  a  strain  of  0.1.  The  cartridge  specimen  loss  coef¬ 
ficient  remains  constant  and  the  compression  specimen  loss  coefficient 
falls  above  a  strain  value  of  0.02 

Based  on  these  results,  the  generally  best  operating  condition  is  wil'.h 
a  double  amplitude  strain  of  0.02.  Of  course,  the  question  of  surviv¬ 
ability  and  life  must  be  assessed  further  before  any  overall  conclu¬ 
sion  about  operating  conditions  are  drawn.  For  shear  elements  the 
temperature  developed  due  to  self-heating  are  unlikely  to  be  a  problem, 
but  fatigue  life  for  significant  strain  amplitude  at  frequencies  in 
the  range  100  to  1000  Hz  must  be  determined  to  be  acceptable.  For 
compression  elements  and  probably  for  cartridge  elements,  the  temper¬ 
atures  developed  are  likely  to  be  a  problem  in  addition  to  fatigue 
life. 

It  is  of  particular  interest  that  the  two  different  cartridge  geom¬ 
etries  show  such  similarity  of  loss  coefficient  dependence  on  strain. 

It  is  also  of  importance  to  note  that  the  cartridge  geometry  indicates 
the  least  overall  sensitivity  of  loss  coefficient  to  strain,  and 
achieves  a  loss  coefficient  value  of  0.2  or  greater  for  all  strains 
above  0.003. 

High  Temperature  Results 

As  expected  the  test  results  show  a  significant  influence  of  tempera¬ 
ture  on  stiffness  and  damping  and  on  the  material  properties  inferred 
from  the  shear  specimen  results.  The  method  of  reduced  variables  (the 
WLF  law)  was  discussed  in  a  previous  section  of  this  report  as  a  means 
of  unifying  the  effects  of  temperature  and  frequency. 

The  availability  of  test  results  for  material  properties  over  a  range 
of  temperatures  now  allows  the  method  of  reduced  variables  to  be  tested 
for  this  material. 

As  identified  the  method  of  reduced  variables  requires  knowledge  of 
the  characteristic  temperature.  No  information  is  available  as  to  the 
characteristic  temperature  or  transition  temperature  for  the  material 
tested  and  it  has  been  necessary  to  seek  to  infer  the  characteristic 
temperature  from  the  test  data.  To  do  this  the  method  has  been  applied 
for  a  range  of  characteristic  temperatures.  Reduced  storage  and  loss 
moduli,  referred  to  32  C,  have  been  calculated  for  each  test  data  point 


for  each  temperature  tested.  The  result  is  a  set  of  discrete,  reduced 
storage  and  loss  modulus  values,  covering  a  wide  range  of  reduced  fre¬ 
quency  values  (from  11  to  1600  Hz),  referred  to  32  C.  Power  law  re¬ 
lationships  between  reduced  moduli  and  reduced  frequency  have  been 
fitted  to  the  data  and,  as  a  measure  of  the  consistency  of  the  reduced 
data,  the  rms  deviation  between  test  data  and  the  power  law  lines  have 
been  calculated.  The  characteristic  temperature  has  then  been  sought 
as  that  value  which  minimizes  the  rms  deviation. 

In  Figures  116  and  117  rms  deviation,  for  reduced  storage  and  loss 
moduli,  respectively,  arc  plotted  against  characteristic  temperature. 

The  deviation  for  storage  modulus  shows  a  minimum  at  a  characteristic 
temperature  oi  -50  C.  Tlie  deviation  for  loss  modulus  however  shows 
no  sign  of  possessing  a  minimum.  There  is  a  continual  decrease  in 
deviation  as  the  characteristic  temperature  is  increased  from  -100  to 
+32  C. 

Accepting  the  storage  modulus  minimum  as  indicating  the  transition 
temperature  and  the  loss  modulus  data  as  anomalous,  reduced  moduli 
have  been  plotted  against  reduced  frequency  for  a  characteristic  tem¬ 
perature  of  -50  C  in  Figures  118  and  119.  These  plots  reflect  the 
minimization  results.  The  reduced  storage  modulus  follows,  with 
reasonable  consistency  a  single  straight  line  on  a  log-log  plot;  the 
reduced  loss  modulus  data  remains  a  series  of  individual  curves 
corresponding  to  the  four  different  temperatures  tested. 

What  are  the  explanations?  In  summary,  either  polybutadiene  does  not 
fall  into  the  category  of  materials  satisfying  the  reduced  variables  law 
or  some  other  parameter  bt'sides  frequency  and  temperature  is  affect¬ 
ing  the  damping  and  loss  modulus  data.  In  the  light  of  the  preceding 
discussion  of  high  dissipation  results  the  possibility  exists  that 
the  influence  of  strain  softening  is  reflected  even  in  the  low 
(reference)  dissipation  results  for  the  shear  specimens. 

Since  increasing  frequency  at  constant  dissipation  means  a  reduction 
in  strain,  the  low  dissipation  test  data  for  damping  is  likely  to 
reflect  the  effects  of  strain  if  these  are  significant.  If  now  we 
review  the  effects  of  strain  on  damping  for  shear  specimens  (Fig.  Ill), 
it  is  apparent  that  decreasing  strain  will  also  decrease  damping. 

Thus,  on  the  basis  of  strain  effects  alone,  "low  dissipation"  damping 
would  be  expected  to  fall  with  increasing  frequency.  This  anticipated 
trend  is  a  plausible  reason  for  the  inconsistency  of  loss  modulus  data 
and  the  method  of  reduced  variables  . 

Stiffness  is  seen  to  increase  with  decreasing  strain  and,  therefore, 
with  increasing  frequency  at  constant  dissipation.  This  is  a  trend’ 
which  is  consistent  with  the  method  of  reduced  variables,  one  of  whose 
basic  tenets  is  that  an  increase  in  frequency  has  an  effect  correspond¬ 
ing  to  a  decrease  in  temperature. 


d 


d 


m 

•H 

0) 

u 

o  ^ 


M 

u 


CO  iH 
f>  W 

CO  u 

iH  60 

0  -H 
'O  w 
o  ^ 

®  c 

CO  d) 
CO  g 
O  -H 
i-j  a 
0) 

d  P. 

•H  C/D 

P  H 

o  cd 

•rl  0) 
Cd  C/D 


o 


60 

•H 


180 


r 


L 

k 

If 

iw 

iv 

L 

m 

jt 


U) 

> 


CO 

;;3 

rH 

0 

•o 

I 

0) 

faC 

ctf 

M 

O 

u 

w 


TJ 

0) 

o 

0 

OJ 

(U 

pc; 


0) 


<y  u 
o 

d  o 

0)  in 


CS3 
•H 
U 

J> 

0) 
o 
0 

4> 
pt; 

'S-a 

TS  0) 

o  a 

^  0 

S'S 

s  »i 


I 


>> 

a 

S  _ 

d  <u 
cr  H 

0) 

M  o 

iH 


00 


Pd 


k 


litt 

if 


181 


Temperature,  32  C.  Char- 
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acteristic  Temperature,  -50  C 


Before  any  conclusion  can  be  drawn  with  respect  to  the  validity  of  the 
reduced  variables  method  for  the  particular  material  it  is  necessary 
to  obtain  further  shear  specimen  values  at  amplitudes  sufficiently  low 
to  exhibit  strain  independence  of  the  material  properties.  In  addition 
data  at  constant  amplitude  should  be  obtained  to  eliminate  this  influence 
of  strain. 

With  the  present  material  property  data  for  polybutadiene  there  is 
limited  benefit  in  applying  the  reduced  variables  method  as  a  method 
of  combining  the  effects  of  frequency  and  temperature.  Since  data  is 
available  at  a  range  of  temperature  values  of  interest  it  has  been 
summarized  and  stored  directly  in  that  form. 

For  application  the  curves  of  Figure  120  have  proved  most  usefui.  In 
these  curves  storage  and  loss  modulus  at  100  and  1000  Hz  are  plotter 
against  temperature.  Values  at  these  two  frequencies  may  therefore  be 
read  for  any  temperature  between  32  and  93  C .  To  calculate  stiffness 
and  damping  at  these  two  frequencies  for  a  particular  element  the 
appropriate  geometric  relationships  are  applied,  and  plots  of  stiff¬ 
ness  and  damping  against  frequency  are  made  which  exhibit  a  straight 
line  variation  on  log- log  paper.  This  is  how  the  predicted  data  for 
cartridge  and  cylindrical  specimens  presented  in  the  preceding  Section 
were  generated. 

Consistency  of  Phase  II  And  Phase  III  Results 

To  provide  a  check  on  consistency,  the  material  property  data  obtained 
from  Phase  II  tests  (as  described  earlier  in  this  report)  has  been 
compared  with  the  Phase  III  material  property  values  of  Table  X  for 
32  C.  The  comparison,  presented  in  Table  XIV  shows  that  there  is  good 
consistency  between  the  two  test  series  for  storage  modulus,  but  that 
the  loss  modulus  values  obtained  from  the  most  recent  series  of  tests 
are  substantially  higher  than  the  Phase  II  values. 


TABLE  XIV.  COMPARISON  OF  PHASE  II  AND  PHASE  III  VALUES 
FOR  STORAGE  AND  LOSS  MODULUS 


Test 

Series 

g' 

100  Hz 
N/m^ 

G' 

1000  Hz 
N/m^ 

g" 

100  Hz 
N/m^ 

G 

1000  Hz 
N/m^ 

Phase  II 

1.353x10^ 

2.067x10^ 

2.021x10^ 

1.684x10^ 

Phase  III 

1.369x10^ 

2.189x10^ 

3.659x10^ 

2.728x10^ 

The  reason  for  the  discrepancy  in  loss  modulus  is  the  simplified,  more 
accurate  test  method  for  material  properties  employed  in  the  most  recent 
series  of  tests.  By  eliminating  the  preload  piston  from  the  shear  tests 
a  source  of  error  has  been  eliminated  -  an  error  which  had  been  exaggerated 
(rather  than  reduced)  by  the  "corrections"  applied  to  account  for  it  in 
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reducing  the  Phase  II  data*.  Contrary  to  expectations  based  on  an  anal¬ 
ysis  of  the  mechanical  system,  assuming  rigid  system  elements  (ref.  2), 
the  effect  of  the  piston's  presence  is  to  reduce  the  apparent  damping 
as  described  earlier  in  this  report.  Since  the  correction  was  based 
on  the  analysis  of  the  mechanical  system,  the  correction  was  applied 
with  the  wrong  sign.  The  reasons  for  this  anomalous  effect  of  the 
piston  are  not  clearly  understood,  but  may  well  result  from  flexibility 
of  such  elements  as  the  piston  itself.  The  fact  that  the  present  test 
data  for  shear  specimens  is  more  reliable  is  not,  however,  in  doubt. 


Comparison  of  Prediction  And  Test  Data 

The  test  data  for  cartridge  stiffness  and  damping  and  for  compression 
element  stiffness  and  damping  at  different  temperatures  and  the  com¬ 
parison  with  prediction  methods,  has  provided  a  basis  for  assessing 
the  worth,  in  particular  applications,  of  the  prediction  methods  dis¬ 
cussed  earlier 

For  cartridge  stiffness  and  damping  predictions,  the  available  methods 
are  based  on  static  stiffness  analysis.  Two  simple  methods,  the  Beam- 
Column  analysis  and  the  method  of  Gobel  predict  stiffness  and  damping 
which  generally  bracket  the  test  data. 

The  difference  in  results  of  these  two  prediction  methods  is  close  to 
a  factor  of  two.  The  limitation  in  state-of-the-art  prediction  for 
this  geometry  may,  therefore,  be  summarized  as  follows: 

Using  static  analysis-based  predictions,  dynamic  performances  of  a  ring 
cartridge  may  be  predicted  with  a  band  of  uncertainty  which  represents 
a  factor  of  two  in  stiffness  and  damping;  if  this  geometry  is  to  be 
employed  in  a  particular  application,  it  should  be  ensured  that  satis¬ 
factory  performance  can  be  obtained  for  values  of  stiffness  and  damping 
which  cover  this  band  of  uncertainty  in  cartridge  predictability. 

For  cylindrical  compression  specimens  there  exists  a  static  prediction 
method,  employing  a  constant  shape  factor  coefficient  value  (g)  of  two. 
Based  upon  this  method,  dynamic  stiffness  and  damping  may  be  predicted 
using  the  frequency  dependent  material  properties  of  storage  and  loss 
modulus.  Alternatively,  using  previously  generated  dynamic  test  data 
for  cylindrical  cartridge  specimens  with  various  d/h  ratios,  a  function 
has  been  developed  empirically  which  provides  frequency  dependent  values 
for  shape  factor  coefficient.  These  frequency  dependent  coefficients 
may  again  be  used,  in  conjunction  with  the  frequency  dependent  storage 
and  loss  module,  to  predict  stiffness  and  damping. 

For  the  Phase  II  test  data,  the  static  based  method  generates  dis¬ 
crepancies  of  as  much  as  50  percent  relative  to  the  test  data. 


For  the  Phase  III  test  data,  at  temperatures  of  66  and  93  C,  the  two 
prediction  methods  produce  discrepancies  of  similar  magnitude.  Figures 
63  through  66  show  the  discrepancies  graphically.  The  worst  discrep¬ 
ancies  between  the  prediction  methods  and  test  data  points  for  all 
test  condition  of  preload  and  temperature  are  summarized  in  Table  XV, 
For  stiffness  the  worst  discripancies  vary  from  8  to  25  percent;  for 
damping  the  discrepancies  are  in  the  range  20  to  30  percent,  with 
measurements  lower  than  predicted. 


TABLE  XV.  CYLINDRICAL  COMPRESSION  SPECIMEN 

[Maximum  Discrepancies  Between  Measured  And 
Predicted  Stiffness  And  Damping  Value] 


Fig , 
No. 

Temp  , 

(C) 

Pre¬ 

load 

(7o) 

Stiffness 

Damping; 

Discrepancy 

Percent^' 

Discrepancy 

Percent 

“  3'  =  2 

3'  =  f(cD) 

-- — — j. — 

3=2 

T'  =  f(a)) 

63 

66 

2.5 

-8  .8** 

+  10.8 

-23.7 

-25  .8 

64 

66 

5 

+  12 ,1 

+13.4 

-28.6 

-30.8 

65 

93 

2.5 

+13,0 

+14.4 

-20.7 

-23.5 

66 

93 

5 

+22.9 

+25.0 

-21.8 

-25,2 

*Percentage  referred  to  predicted  value. 

’'^''Negative  sign  indicates  measured  data  lower  than  predicted. 

_  Positive  sign  indicates  measured  data  higher  than  predicted. 


In  terms  of  performance,  there  is  nothing  to  choose  between  the  two 
approaches  for  cylindrical  compression  elements.  When  either  method 
IS  to  be  used  in  design  of  a  system  component,  it  is  necessary  to 
ensure  that  satisfactory  system  performance  can  be  achieved  allowing 
ror  uncertainties  of  the  magnitude  indicated  in  Table  XV. 


Test  Method 

The  ability  of  the  base  excitation  resonant  mass  method  to  provide 
stiffness  and  damping  data  for  a  variety  of  geometries  and  over  a 
range  of  temperatures  and  dissipation  levels  has  been  convincingly 
demonstrated  by  the  results  presented  in  this  report.  It  has  bLn 
determined  that  the  "cleanest"  test  condition  occurs  when  the  resonant 
mass  IS  supported  solely  by  the  test  specimen,  without  either  of  the 
lower  or  upper  loading  and  unloading  pistons.  This  simple  test  con¬ 
dition  introduces  no  uncertainty  with  respect  to  tare  stiffness  and 
damping . 


The  simple  test  condition  with  no  pistons  can  be  achieved  for  the 
cartridge  and  shear  specimens,  without  danger  of  putting  the  elastomer 
metal  bond  in  tension.  The  compression  specimen  as  presently  designed 
requires  the  preloading  effect  of  the  piston  to  maintain  the  bond  in 
compression.  A  modified  design  can  be  envisaged  in  which  a  self¬ 
preloading  mechanism  is  employed. 

Future  Directions 

The  results  presented  in  this  report  make  clear  the  need  for  further 
work  to  improve  fundamental  understanding  and  applications  capability. 

The  influence  of  strain  amplitude  as  a  parameter  independent  of  theruu  ’ 
dissipation  level  must  be  further  investigated  as  a  function  of  geometry 
for  polybutadiene  and  for  different  materials.  The  present  results  in¬ 
dicate  that,  over  a  significant  range  of  strains,  stiffness  and  damping 
are  governed  by  strain  rather  than  frequency.  This  must  be  further 
tested.  In  addition,  the  conditions  of  strain  for  which  the  stiffness 
and  damping  are  independent  of  strain  must  be  investigated,  and  related 
to  geometry.  It  is  anticipated  that  tests  at  very  low  amplitudes  as 
well  as  high  amplitudes  will  be  required.  The  means  to  present  the 
resultant  information  in  a  form  suitable  for  use  in  design  must  also 
be  developed.  The  information  will  be  largely  empirical,  until  the 
strain  softening  effect  becomes  better  understood,  and  test  methods 
must  be  evolved  which  can  rapidly  define  the  effect  for  different 
materials  and  geometries. 

To  improve  understanding  of  the  importance  of  thermal  dissipation  and 
its  relationship  to  strain  softening,  thermal  properties  of  conductivity 
and  specific  heat  must  be  determined  for  the  material  being  tested  and 
included  correctly  in  the  thermoviscoelastic  model.  Direct  measurements 
of  temperature  within  the  elastomer  material,  and  comparison  of  these 
measurements  with  predictions  will  also  provide  improved  understanding 
of  thermal  dissipation  effects. 

When  strain  softening  is  better  understood,  and  stress  independent 
results  are  available,  the  method  of  reduced  variables  as  a  means  of 
conveniently  presenting  elastomer  properties  must  be  further  evaluated. 

As  comparisons  show,  the  uncertainty  in  prediction  ability  for  a  car¬ 
tridge  and  of  cylindrical  compression  elements  when  static  geometry 
effects  are  applied,  while  significant,  may  be  acceptable  for  some 
applications.  The  degree  of  predictive  uncertainty  should  be  in¬ 
vestigated  further  for  other  materials.  The  need  for,  and  potential 
benefits  of,  improved  viscoelastic  analysis  methods  must  be  identified 
and,  if  the  need  is  shown  to  exist,  appropriate  methods  of  viscoelastic 
continuum  analysis  should  be  developed. 

Tests  over  similar  parameter  ranges  for  geometries  similar  to  those 
presently  reported  should  be  performed  using  other  materials,  including 
materials  selected  for  known  differences  in  hardness  and  loss  coef¬ 
ficient.  These  tests  will  provide  the  basis  for  generalization  of^  or 
exception  tc^  the  results  for  polybutadiene. 


In  the  case  of  a  cartridge,  the  effect  of  a  rotating  load  as  opposed 
to  a  uniaxial  load  must  be  investigated.  The  anticipated  application 
of  the  cartridge  element  is  in  rotating  machinery  where  the  load 
vector  will  be  rotating.  It  is  important  to  determine,  under  con¬ 
trolled  conditions,  whether  there  is  a  significant  difference  in  the 
stiffness  and  damping  under  the  action  of  rotating  and  uniaxial  dynamic 
loads,  and  to  establish  the  degree  of  uncertainty  in  prediction  for  a 
rotating  load  environment. 

An  important  geometry  to  be  tested  is  the  0-ring,  This  is  a  compact, 
commerc ia 1 ly-ava i lab le ,  simple,  elastomer  element.  0-rings  are  already 
in  use  as  seals  and  flexible  mounts  in  squeeze-film  dampers  where  their 
stiffness  properties  in  particular  are  important.  Their  use  as  flexi¬ 
ble  mounts  independent  of  squeeze  film  dampers  is  anticipated,  in 
which  case  their  damping  properties  will  also  be  important. 

When  testing  0-rings,  and  other  geometries,  it  is  important  to  preserve 
the  approach  which  has  been  inherent  in  the  tests  and  analysis  up  to 
this  point;  that,  for  improved  understanding,  tests  in  shear,  for 
general  material  properties,  must  be  performed  on  the  same  material  as 
is  employed  in  the  geometry  of  interest.  In  this  way,  meaningful  pre¬ 
dictive  methods  which  relate  stiffness  and  damping  for  storage  and  loss 
moduli  can  be  evolved. 

Analytical  studies  for  representative  machinery  classes  should  be 
performed  to  establish  the  needs  for  damping  elements  and  the  ability 
and  limitations  of  elastomers  in  filling  these  needs. 

Application  tests  should  be  performed  to  demonstrate  the  vibration 
control  capabilities  of  elastomer  elements.  Following  successful 
application  tests,  the  questions  of  life  and  survivability  should  be 
investigated . 
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appendix  a 


THERMAL  DISSIPATION  IN  AN  ELASTOMER  ELEMENT  UNDER  UNIAXIAL  STRESS 

As  a  result  of  imposed  conditions  an  elastomer  element  in  f 
o  a  right  circular  cylinder,  is  undergoing  time  varying  changes  in'" 
temperature,  internal  stress  and  displacement.  The  cylfnSorfras 
shown  in  Figure  121.  It  is  assumed  that  stresses  ar^  unia^iaK 

At  any  point 


a 

z 

U 

z 

=  U(z,t)e^“’' 

(73) 

T 

z 

=  Tjz.t) 

(74) 

where  is  the  normal  axial  stress 

U^  is  the  axial  displacement;  is  temperature; 

PKpre.slons  (72)  and 

a^and  U  are  complex  nnmbera  containing  both  amplitude  and  phase  informa- 

It  is  assumed  that  the  amplitude  functions  a,U  are  undergoing  much 
slower  changes  with  time  than  a  and  U  b  e  uiuen 

z  z  * 


Axial  equilibrium  requires: 

^  °  (75) 

where  Z  is  the  axial  coordinate  and  the  inertia  of  the  elastomer 
material  is  neglected,  «oi.umer 

The  stress  amplitude  o  is  evaluated  as: 


a  -  3G*(Z,a))  ^  (75) 

and^c'^  c"  ^  ^  ^  (‘^>T);  T  is  the  cycle  average  temperature; 

and  G  ,  G  are  the  storage  and  loss  moduli.  T  is  a  function  of  Z. 

The  instantaneous  energy  equation  is: 
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P  c  ~  =  K  +  S 
P  dZ^ 


.  .dU 

where  S  =  Re  (a  )  Re  :^  ( 

2-^  dt  VdZ 


Now,  assuming  temperature  to  be  changing  slowly  with  time,  we  can  take 
a  cycle  average  as  follows  to  get  T: 


T  = 


=  ~ 

2n  J 


.t+Zn/oi 


T^(t  )dt 


thus  the  energy  equation  becomes: 


p  c  ^  =  K  4  .  f 

p  3t  ^^2  2tt 


S(t')dt' 


Now,  if 


t+(2rr/ao) 

I  =  \  S(t^)dt 


I  '  4 


where  the  superscript  ~  implies  the  complex  conjugate, 


,  ~  -icut\  S  / .  „  icDt 
+  a  e  j  ^  (^loiJ  e 


-ICDt  \  ,  ^ 


-Ztt/o)  - 


ico  I  ^  2i(ot  ^_~dU  -Zicot'L  , 


az  ^  °  az  az  ^  az 


^  Is.  -  .M' 
2  r  dz  azJ 


(continued  on  next  page) 
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■is.  -jp  M  su  _  ^  du 

2  3Z  9Z  az 


3in  ^ 

G  G 

M  su 

3Z  BZ 

,  2 

J 

Thus,  the  final  form  of  the  energy  equation  is: 


^  P  St  2  2  \l  BZ  / 


and  the  equilibrium  equation  is: 

a  r.  , 


T  r*  tv  \  3U  '] 

az  J  - 


It  is  straightforward  to  perform  a  similar  development  for  elements 
undergoing  shear  deformation  (as  shown  in  Fig.  122),  in  this  case  the 
energy  equation  is: 

-  >=  J  "  (M) 

where  U  is  the  amplitude  of  shear  deformation  (normal  to  Z)  and  the 
equilibrium  equation  is: 

^  [g*  (Z.O))  ||]  =  0  (86) 

To  further  streamline  Equation  (83)  through  (86)  to  a  single  pair  of 
equations,  we  can  write  (reorganizing  slightly)  for  shear  or  compres¬ 
sion  elements: 

Energy  equation: 
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Equilibrium  equation; 


(88) 


where  y  1  for  a  shear  element,  3  for  a  compression  element. 


s 
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APPENDIX  B 


DATA  REDUCTION  ANALYSIS 

For  the  Bese  ExCj.i.ation  Resonant  Mass  test  method ,  test  data  consists 
of  accelerometer  and  displacement  sensor  readings.  The  primary 
objective  of  data  reduction  is  to  interpret  these  signals  in  terms  of 
stiffness  and  damping.  The  structure  of  this  data  and  the  method  of 
reduction  can  best  be  defined  in  relation  to  a  simple  spring-mass- 
damper  system  subjected  to  base  excitation,  as  shown  In  Figure  123. 

The  equation  of  motion  for  the  mass  is; 


where  is  displacement  of  mass 

is  displacement  of  base. 


We  will  assume  that  the  motions  of  base  and  mass  are  harmonic;  so  that 


Y  =  Real  {Y*e^“^} 

m  L  m  j 

Y  -  Real  {y* 

o  L  o  J 


and  the  equation  of  motion  may  be  rewritten: 

-Jm  +  ica  C(Y*-Y*)ei»*  -  0  (92) 


Dividing  through  by  we  obtain 


2  * 
CO  M  Y 


O'!  +  iK2>  = 


where  =  co  C 


Now  if  cp  is  the  phase  angle  between  and  ^we  can  write 


Y*|cos  cp  —  i  Y*  sin  cp 

m  111 
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OUTPUT 


j/ 

_J  INPUT. 


Model  of  Base  Kxcitation  Resonant 
Mass  System  for  Data  Redncticni 
Analysis 


(96) 


(K  +  iK  )  =  -  cos  cp  -  i  Of  sin  co: 

1  2  O'  cos  cp“icv  sincp~l 

|Y*| 


where  a  is  the  ti.ansmissibility-^-y'-:  Thus 


=  cp  M  (a  cos  (p  i  g  sin  cpXcy  cos  cp  “  1  +  i  a  sin 

1  2  2  2  2 

(o'  cos  cp  1)  +  O’  sin  cp 


CD  M  [cy^  —  a  cos  cp] 


2a  cos  cp  +  1 


tr  _  _..P  M  a  sin  cp  _ 

-  (99) 

a  ~~  2a  cos  cp  +  1 

Thus,  if  we  determine  a,  the  transmissibility,  or  ratio  of  output 
amplitude  |Yjjjj  to  input  amplitude  |Y*  |  ,  and  cp  the  phase  angle  between 
input  and  output  amplitudes,  we  can  directly  evaluate  stiffness  and 
damping  (Kj^  and  K2). 

these  expressions  are  based  on  the  simplest  one"degree“Of" 
freedom  model  with  no  correction  for  tare  damping  or  stiffness.  This 
is  in  contrast  to  the  expressions  used  in  Reference  2,  which  made 
corrections  for  effects  of  damping  in  preload  and  load  support  pistons. 
The  simplicity  of  the  present  treatment  results  from  the  simplicity  and 
accuracy  of  the  test  procedure  performed  for  the  majority  of  the 
tests,  in  which  no  external  restraint  of  any  kind  was  applied  to  the 
supporting  mass.  Only  for  the  compression  specimens  was  the  preload 
piston  employed.  These  expressions  are,  therefore,  fully  accurate  for 
shear  and  cartridge  specimens  and  subject  to  some  uncertainty  only 
in  the  case  of  the  compression  specimens. 
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NOMENCLATURE-: 


A  stressed  area,  ni“ 

B 

^^^2  coefficients  of  power  law  relationship  K  =  A.oo 

^2  2 
=  A^ci'j  ,  N/m 

"^o ’"^i’^2 ’^^3  coefficients  of  force  in  generalized  force-displace¬ 
ment  relationship 

^0  factor  by  which  frequency  is  multiplied  to  get  reduced 

frequency  at  temperature  0  ,  corresponding  to 
temperature  6  ^ 

2 

B  bulk  modulus,  N/m 

b  width  of  shear  element  in  direction  perpendicular  to 

direction  of  loading,  m 

B 

exponents  in  power  law  relationship  =  A  O)  K  = 

B  112 

A^co  ^ 

^o»^1»^2’^3  coefficients  of  displacement  and  its  time  derivative 
in  generalized  force-displacement  relationship 

specific  heat,  J/kg  C 

damping  coefficient  of  damping  element  in  triparameter 
model,  Ns/m 

diameter  of  cylindrical  compression  element,  m 
elemental  area,  m^ 

elemental  forces  in  Beam-Column  analysis,  parallel  and 
perpendicular  to  Beam-Column  axis,  respectively,  N 

inner  and  outer  diameters  of  cartridge  element,  m 

energy  loss  per  cycle  due  to  hysteresis,  J 

elemental  time  in  cycle  average  integral 

elemental  angle 

Young's  modulus,  N/m^ 

effective  Young* s  modulus  for  cartridge  element,  N/m^ 


T1 

D 

da 


df  ,df 
c  s 


D.  ,D 
1  o 


s 

dt ' 

de 

E 
E 


eff 


G  ,G 


K 


K 


Tl’  T2 
K, 


K. 


1,1. 


force  exerted  by  elastomer,  N 
frequency,  Hz 

maximum  force  in  hysteresis  loop,  N 

factor  used  by  Gobel  to  correct  for  shape  effects  in 
cartridge  elements  of  finite  length 

reduced  frequency,  Hz 

complex  force  amplitude,  N 

O 

shear  modulus,  N/m 
complex  modulus  G  ^  4-  iCf  N/m 

2 

storage,  loss  moduli  in  shear,  N/m 
effective  storage  modulus  in  shear,  N/m^ 
effective  loss  modulus  in  shear,  N/m^ 

n 

Static  shear  modulus,  N/m 
reduced  shear  modulus,  N/m^ 

height  of  compression  element;  thickness  of  shear 
element  (strained  dimension  for  compression  and  shear 
elements) 

(Ixj/co)  X  cycle  average  of  dissipation  function,  S, 
J/m^/cycle 

thermal  conductivity,  watts /cm  C 
complex  stiffness,  +  iK^ ,  N/m 
radial  stiffness  of  cartridge,  N/m 
stiffness  values  in  triparameter  model,  N/m 
stiffness,  N/m 
damping,  N/m 

length  of  cartridge  element,  m 
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M 

n,N 

r" 


supported  mass,  kg 

number  of  elements  in  test  specimen 

mean  radius 

right-hand  side  of  energy  equation  in  numerical 
solution  at  time,  t,  J/m^S 

inner  and  outer  radius  of  cartridge,  m 

dissipation  function,  j/m^S 

shape  factor  -  ratio  of  one  loaded  area  to  free  area 
temperature,  C 
t ime ,  s 

dummy  time  for  cycle  average  integral,  s 

ambient  temperature  in  metal  to  which  elastomer  is 
bonded,  C 

instantaneous  temperature  at  axial  coordinate  Z,  C 
axial  displacement,  m 

energy  stored  in  element  at  maximum  displacement,  J 

axial  displacement  (instantaneous),  m 

length  of  shear  element  in  direction  parallel  to 
loading,  m 

deflection  of  elastomer  element,  m 

complex  amplitude  of  elastomer  deflection,  m 

absolute  displacement  of  supported  mass,  m 

complex  displacement  amplitude  of  supported  mass,  m 

absolute  displacement  of  base,  m 

complex  displacement  amplitude  of  base,  m 

maximum  deflection  of  elastomer  element,  m 

axial  coordinate,  m 


'■% 

J 

J 

i 

J 

J 

i 

J 

i 

i 


mi 


d 
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transmissibility,  ratio  of  amplitude  magnitudes:  output 
acceleration  to  input  acceleration 

shape  factor  coefficient 

shape  factor  coefficient  for  use  in  calculating 
stiffness 

shape  factor  coefficient  for  use  in  calculating 
damping 

factor  equal  to  1  for  shear  elements,  3  for  conorassion 
elements 

time  step  in  numerical  solution  of  energy  equation,  c 

loss  coefficient, 

angle  in  cartridge  analysis 

absolute  temperature  in  reduced  variables  analysis,  K 

characteristic  temperature  for  elastomer  material,  K 

transition  temperature  for  elastomer 

density  of  elastomer,  kg/m^ 

complex  stress  amplitude,  N/m^ 

normal  axial  stress  (instantaneous),  N/m^ 

phase  angle  between  motions  of  base  and  supported  mass 

frequency,  rad/sec 
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